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Abstract 

We study the complexity of a class of problems involving satisfying constraints which remain the 
same under translations in one or more spatial directions. In this paper, we show hardness of a classical 
tiling problem on an iV x iV 2-dimensional grid and a quantum problem involving finding the ground state 
energy of a 1-dimensional quantum system of N particles. In both cases, the only input is N, provided 
in binary. We show that the classical problem is NEXP-complete and the quantum problem is QMAexp- 
complete. Thus, an algorithm for these problems which runs in time polynomial in N (exponential in 
the input size) would imply that EXP = NEXP or BQEXP = QMAexp, respectively. Although tihng 
in general is already known to be NEXP-complete, to our knowledge, all previous reductions require 
that either the set of tiles and their constraints or some varying boundary conditions be given as part of 
the input. In the problem considered here, these are fixed, constant-sized parameters of the problem. 
Instead, the problem instance is encoded solely in the size of the system. 

1 Introduction 

One perennial difficulty with practical applications of hardness results is that the practically interesting 
instances of a hard language may not themselves form a hard class. One approach to solving this problem 
is the difficult theory of average-case complexity IILev861IBDC G891. in which one can show that "typical" 
cases of some language are hard. In this paper we take a different approach. In many cases, practically 
interesting instances possess some shared property, such as a symmetry, that distinguish them from the 
general instance and might, in principle, make those instances easier. We will study such an example and 
show that, even in a system possessing a great deal of symmetry, it is still possible to prove a hardness result. 

Specifically, we consider the related problems of determining whether there is a possible tiling of an 
r-dimensional grid with some fixed set of classical tiles and of finding the lowest energy state (or ground 
state) of a quantum system involving interactions only between neighboring particles on an r-dimensional 
grid. The ground state energy of a system is considered one of the basic properties of a physical system, 
and over the last few decades, physicists have developed a number of heuristics that have been successful 
in finding the ground state energy in many special cases. On the other hand, in earlier work IIAGIK07II . we 
have shown that in the most general case, even in a 1-dimensional quantum system, finding the ground state 
is a computationally difficult problem (modulo the usual complexity-theoretic assumptions). However, the 
construction presented in [ AGIK071 involves a system which is completely unnatural from a physical point 
of view. The most interesting physical systems frequently possess an additional symmetry: translational 
invariance. In this paper, we will show that even a 1-dimensional translationally-invariant system can be 
hard. 

One interesting feature of our proof which may have more general applicability is that the only free 
parameter for the language we consider is the size of the system. This is frequently the case for interesting 
systems: there is a basic set of rules of constant size, and we wish to study the effect of those rules when the 
system to which the rules apply becomes large. In practice, many such systems seem difficult to solve, but 
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it is hard to see how to prove a complexity-theoretic hardness result, since that requires reducing a general 
problem in some complexity class to the language under consideration, and there doesn't seem to be room 
in the language to fit all the needed instances. Usually, this difficulty is circumvented by modifying the 
problem slightly, to add additional parameters in which we can encode the description of the instance we 
wish to simulate. 

To illustrate, let us present the classical tiling problem we study in this paper: We are given a set of 
square tiles which come in a variety of colors. The area to be tiled is a square area whose size is an integer 
multiple of the length of a tile. We are given horizontal constraints indicating which pairs of colors can be 
placed next to each other in the horizontal direction and another set of constraints in the vertical direction. 
We specify a particular color which must go in the four comers of the grid. The description of the tile colors, 
placement constraints and boundary conditions are fixed for all inputs of the problems. The input is just a 
number N written in binary and we wish to know whether an x grid can be properly tiled given these 
constraints. We show that this problem is NEXP-complete. Note that the input in this case is size log N, so 
an algorithm to solve our tiling problem that runs in time polynomial in would imply that NEXP = EXP. 
While it is possible that P 7^ NP and yet NEXP = EXP, this seems unlikely to be the case. 

This version of tiling is equivalent to the more common Wang Tiles | Wan601 in that any set of tiles can be 
transformed into a set of Wang Tiles (and vice versa) such that there is a one-to-one correspondence between 
valid tilings on an x grid. For an infinite grid, the problem is undecidable. Intuitively, it makes sense 
that it is also hard (for some sets of tiles) for a finite grid, since there are exponentially many possible tilings, 
and it is impossible to tell locally whether a given partial tiling can be extended indefinitely. Indeed, there 
are prior results showing that related tiling problems are NEXP-complete, but to our knowledge, all previous 
reductions require that either the set of tiles and their constraints or some varying boundary conditions be 
given as part of the input l|LP971IBoa97l| . For instance, one may specify the placement of some number of 
tiles and ask whether it is possible to extend that partial tiling to a tiling of the full square. Even though that 
problem had been proven hard, the more natural problem of whether it is possible to efficiently find a tiling 
of the empty grid remained open. 

The basic idea of our construction is to reduce from an instance x of some language in NEXP by 
encoding x in the binary expansion of A^, the size of the grid. It is well known that a finite set of tiling rules 
can be used to implement a universal Turing Machine. We need some way to express the program for the 
Turing Machine to run, and that program must grow with the size of x. Previous constructions managed 
this by resorting to either polylog A^ different tile types or varying boundary conditions to encode x, but 
those are both fixed, constant-sized parameters in our version of the problem. Instead, we use the tiles to 
implement a binary counter which converts A^ into binary and then uses it as an input to a universal Turing 
Machine. 

The other problem we consider is finding the ground state energy of a quantum system. The state of a 
quantum system with A^ qubits is a vector in a Hilbert space of dimension 2^. We will be considering a 
slightly more general version in which an individual particle has its state in a space of dimension d, in which 
case the state of a system of A^ such particles is a vector in a -dimensional Hilbert space. One of the 
postulates of quantum mechanics states that any physical property of a system that can be measured (e.g. 
location, momentum, energy) corresponds to a linear operator. For an A^-particle system, it can be expressed 
as a X matrix over the complex numbers. If the property is measured, then the outcome must be an 
eigenvalue of the corresponding linear operator and the state of the system after the measurement is in the 
eigenspace corresponding to the outcome. Thus, the problem of finding the energy for the lowest energy 
state is the same as determining the lowest eigenvalue for the energy operator (also called the Hamiltonian 
for the system). The difficulty, of course, is that the Hamiltonian matrix is exponentially large in the size A'^ 
of the system. 

We are typically interested in systems whose Hamiltonians are local in that they can be expressed as 
a sum of terms each of which acts non-trivially only on a constant-sized subset of the particles in the sys- 
tem. Although the term "local" does not imply anything about the physical location of the particles, it is 
motivated by the idea that particles only interact when they are physically close to each other. We are there- 
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fore interested in extending this even further and examining particles that are located in a geometrically 
r-dimensional space where only particles within a fixed distance can interact. A particularly natural model 
to consider, then, is a system of particles on an r-dimensional grid, where the terms of the Hamiltonian 
operate only on neighboring pairs of particles in the grid. Note that although the full matrix representation 
of a Hamiltonian is exponentially large in the size of the system, a local Hamiltonian has a compact repre- 
sentation: each term can be expressed as a constant-sized matrix, and there can only be polynomially many 
such terms. 

Kitaev introduced the class QMA, the quantum analog of NP, and showed that the problem of deter- 
mining the ground state energy of a system defined by a local Hamiltonian is QMA-hard IIKSV02II . Thus, 
we do not hope to solve it even on a quantum computer. With an additional promise, the problem is QMA- 
complete: there exist two values a > b, such that a — b > l/poly(A^), where it is guaranteed that the 
ground state energy is at most b or at least a, and one wants to determine only which of the two alternatives 
holds. The problem is still hard even for two-dimensional systems on qubits or one-dimensional systems of 
particles of constant Hilbert space dimension IIOT051IAGIK07II . 

Despite these worst-case results, numerical methods have been successful at determining ground state 
energies for many quantum systems, especially in one dimension. What are the differences between these 
hard QMA-complete problems and the more tractable systems studied by numerical physicists? One fea- 
ture of the QMA-completeness constructions is that the individual terms of the Hamiltonian are position- 
dependent. Essentially, the computation performed by a quantum verifier circuit is encoded into the Hamil- 
tonian so that a low energy state exists if and only if there is a quantum witness that causes a verifier to 
accept. Thus, the terms of the Hamiltonian encode, among other things, individual gates in a quantum cir- 
cuit. In contrast, many quantum systems of physical interest are much more uniform in that they consist 
of a single Hamiltonian term that is simultaneously applied to each pair of neighboring particles along a 
particular dimension. Such a system is called translationally invariant. 

Since highly symmetric systems are rather natural, a number of researchers have studied the com- 
putational power of translationally invariant quantum systems. For instance, UNWOSII gives a 20-state 
translation-invariant modification of the construction from IIAGIK07II (improving on a 56-state construc- 
tion by [JWZOTjI) that can be used for universal 1-dimensional adiabatic computation. These modifications 
require that the system be initialized to a particular configuration in which each particle is in a state that 
encodes some additional information. The terms of the Hamiltonian, although identical, act differently on 
different particles depending on their state. The ground state is therefore degenerate and one determines 
which ground state is reached by ensuring that the system starts in a particular state. Kay pKayOS] gives a 
construction showing that determining the ground state energy of a one dimensional nearest-neighbor Hamil- 
tonian is QMA-complete even with all two-particle terms identical. The construction does, however, require 
position-dependent one-particle terms. Irani has demonstrated ground state complexity in one-dimensional 
translationally-invariant systems by showing that such systems can have ground states with a high degree 
of quantum entanglement ilra09ll . While quantum entanglement is closely related to the performance of 
numerical heuristics in practice, the particular states in this construction are easy to compute. 

In contrast, we show that there exist 1-dimensional translationally-invariant quantum systems with 
nearest-neighbor interactions for which finding the ground state energy is complete for QMAexp> a quantum 
analogue of NEXP. As with the classical result, the only parameter which varies in the language is A^, the 
number of particles, and we must use N to encode the instance from which we wish to reduce. The quantum 
result uses a similar idea to the classical result: we arrange for a control particle to shuttle between the ends 
of the system and count the number of particles. The binary encoding for the number of particles is then 
used as an input to a quantum Turing Machine. 

It is worth noting that the one-dimensional version of the classical tiling problem is very easy: it is in 
P (see section 1531 for the algorithm). That is, it can be solved in a time polylog N, whereas it appears the 
quantum problem takes time exp(A^), even on a quantum computer (unless QMAexp = BQEXP, where 
BQEXP is like BQP, but with exponential circuits). Translational invariance does seem to simplify the 1- 
dimensional classical case, reducing poly(A^) time to polylog(A^) time, but it doesn't help very much in the 
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quantum case. 

Note that the classical tiling problem is a special case of the ground state energy problem for quantum 
systems where the Hamiltonian is diagonal in the standard basis with only 1 or entries. Any ground state 
of such a system is a classical state in which the state of each particle is specified by one of the d possible 
standard basis states, which correspond to the possible tile colors. A pair of tiles {ti,tj) is allowed by the 
tiling rules iff the corresponding \titj) {titj \ term of the Hamiltonian is 0, so that allowed tilings have total 
energy, whereas a forbidden tiling has energy at least 1. 

2 Problems and Results 

Definition 2.1 TILING 

Problem Parameters: A set of tiles T = {ti, . . . ,tm}- A set of horizontal constraints H C T x T such 

that ifti is placed to the left oftj, then it must be the case that {ti, tj) G H. A set of vertical constraints 

V C. T X T such that ifti placed below tj, then it must be the case that (ti, tj) gV. A designated tile ti 

that must be placed in the four corners of the grid. 

Problem Input: Integer N, specified in binary. 

Output: Determine whether there is a valid tiling of an N x N grid. 

Theorem 2.2 TILING is UEXP -complete. 

We give the proof in section [3i The basic idea is that two adjacent comers are used to create a border 
around the perimeter of the grid which allows us to implement special rules at the top and bottom rows. The 
interior of the grid is tiled in two layers, each of which implements the action of a Turing machine. The first 
TM proceeds from top to bottom on layer 1 and the second proceeds from bottom to top on layer 2. The first 
TM takes no input and acts as a binary counter for N steps. The bottom row of the first layer then holds a 
binary number that is Q{N^/''). The rules for the lower boundary are then used to copy the output from the 
binary counter to the bottom row of layer 2, which acts as the input to a generic non-deterministic Turing 
machine. The rules for the top boundary check whether the final configuration on layer 2 is an accepting 
state. 

A common convention for this tiUng problem is to only specify the boundary condition tile in a single 
corner of the grid. This does not work in our case, so we instead use specified tiles in all four corners to mark 
out the boundary of the grid to be tiled. We have considered other versions of the classical translationally- 
invariant tiling problem to understand to what extent the precise definition of the problem is important. The 
boundary conditions are a critical component. We have chosen to fix the tiles used at the 4 corners of the 
square, but there are a number of other possible ways to set the boundary conditions, including periodic 
boundary conditions (so we are actually tiling a torus) and open boundary conditions, where any tile is 
allowed at the edges of the square. The case of periodic boundary conditions is particularly interesting 
because it is truly translationally invariant, unlike our usual formulation where the boundaries break the 
translational symmetry. We show this case is also hard, but with a more complicated reduction than in our 
standard TILING problem. Another variant is to make the problem more similar to the quantum Hamiltonian 
problem by assigning a cost to any pair of adjacent tiles, and allowing the costs to be different from or 1. 
This is like a weighted version of tiling and corresponds to a Hamiltonian which is diagonal in the standard 
basis but does not have any other constraints. 

We have also considered problems with additional symmetry beyond the translational invariance. If we 
have reflection symmetry, then if {ti,tj) G H, then {tj,ti) G H as well, and if {ti,tj) € V, then {ti,tj) € V 
also. That is, the tiling constraints to the left and right are the same, as are the constraints above and below. 
However, if we only have reflection symmetry, there can still be a difference between the horizontal and 
vertical directions. If we have rotation symmetry, we have reflection symmetry and also {ti,tj) € H iff 
{ti,tj) € V. Now the direction does not matter either. These additional symmetries are well motivated from 
a physical point of view since many physical systems exhibit reflection and rotation symmetry. We do not 
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2-D, no symmetry 



2-D, reflection sym. 



2-D, rotation sym. 1-D 



BC on all corners 

unweighted 
weighted 
BC on or 1 corner 
unweighted 
weighted 
Periodic BC 

unweighted 
weighted 



NEXP-complete^ 
NEXP-compIete 



P, uncomputable 
NEXP-complete 



NEXP-compIete 
NEXP-complete 
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P, uncomputable? 

7 



P 



P, uncomputable? 
P, uncomputable? 



P 
P 



P 
P 



P 
P 



P 
P 



P 
P 



P 
P 



Table 1: Summary of the variants of TILING. "BC" is short for "boundary condition." "P, uncomputable" 
means that the associated problem is in P, but an essential parameter of the efficient algorithm we found is 
uncomputable (with a question mark if we are not sure whether it is uncomputable). "NEXP-complete*" 
means complete under an expected poly-time randomized reduction or a deterministic polyspace reduction. 

know how these symmetries affect the complexity of the quantum case, and do not have complete results 
for the classical case. (In particular, we have not solved the weighted variant with reflection symmetry for 
open or four-corners boundary conditions.) Finally, we have studied the one-dimensional version of the 
problem as well as the two-dimensional version. See Table[T]for a summary of our results. Proofs are given 
in Section [51 

Some variants of TILING we consider are easy but in a strange non-constructive sense in that there 
exists A'o G Z"*" U {oo} such that if < A'^o. there exists a valid tiling, and if > A'o, then there does not 
exist a tiling (or sometimes the other way around). However, Nq is uncomputable as a function of (T, H, V). 
These cases are denoted as "P, uncomputable" in Table [T] (with a question mark if we have not been able 
to prove whether A'^o is computable or not). Note that this does not exclude the existence of a (potentially 
slower) algorithm to solve particular instances; indeed, all the classes in Table[T]are included in NEXP. For 
these variants, we know that there is an efficient algorithm, so a hardness result can be ruled out, but since 
the algorithm depends on an uncomputable parameter, it may be that the problem remains hard in practice. 

Now we turn to the quantum problem. First we need to define the class QMAexp- It will be a bit more 
convenient to work with quantum Turing Machines than quantum circuits. The definition is the same as 
QMA except that the witness and the length of the computation for the verifier (which is a quantum Turing 
Machine) can be of size 2" on an input of length n. 

Definition 2.3 A language L is in QMA^xp iff there exists a k and a Quantum Turing Machine M such 
that for each instance x and any on 0(2l^l ) qubits, on input (x, M halts in 0(2l^l ) steps. 
Furthermore, (a) if x G i^yes> 3 such that M accepts (x, with probability at least 2/3. (b) if 
X G Ljio. then V \ip), M accepts {x, {ip)) with probability at most 1/3. 

Definition 2.4 r-DIM TIH (Translationally-Invariant Hamiltonian) 

Problem Parameter: r Hamiltonian terms Hi , . . . , Hr that each operate on two finite dimensional parti- 
cles, specified with a constant number of bits. Two polynomials p and q. 
Problem Input: Integer N, specified in binary. 

Promise: Consider an -dimensional grid of particles and the Hamiltonian resulting from applying Hi 
to each pair of neighboring particles along dimension i. The ground state energy of this system is either at 
mostp{N) or at least p{N) + l/q{N). 

Output: Determine whether the ground state energy of the system is at most p{N) or at least p{N) + 



The following theorem is the main result for the quantum case and shows that the problem will likely 
be, in general, difficult. Note that typically, one is willing to spend time that is polynomial in the size of 



l/q{N). 
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the system (which is in turn exponential in the size of the input). It follows from the result that if there is 
a quantum algorithm that finds the ground state energy in time that is polynomial in the size of the system 
then QMAexp = BQEXP. 

Theorem 2.5 I-DIM TIH is QUk^xp -complete. 

The theorem immediately implies that r-DIM TIH is QMAEXp-complete for any r > 1 since we can 
always take Hi = for i > 2 which results in a system of A^*"^^ independent lines with N particles. We 
prove theorem [23] in section ID If we instead define r-DIM TIH with periodic boundary conditions, we get 
the same result. 

As is common in QMA-completeness results, the construction for Theorem 12.51 creates a Hamiltonian 
whose ground state is a uniform superposition of a sequence of states which represent a particular process. 
A portion of the Hilbert space for the system holds a clock which allows us to control the length of the 
sequence and ensures that the states in the sequence are mutually orthogonal. That is, the t*^ state has the 
form \4)t)\t), where \(j)t) is the t*^ state in the process we wish to simulate, and the overall ground state will 
be \4)t)\t). The size of the system controls the number of time steps for which the clock runs. In the case 
of the construction presented here, the process consists of two main phases. The first phase is the execution 
of a Turing machine which simply increments a binary counter. The clock ensures that this TM is run for 
— 3 steps after which a number that is Q{N^/^) is encoded in binary in the state of the quantum system. 
This state is then used as the input to an arbitrary quantum Turing machine which is executed in the second 
phase. This QTM implements a verifier which is also allowed a quantum witness of length Q{N). Finally, 
there is an energy term which penalizes any non-accepting computation of the verifier. 

As a corollary of Theorem [231 the following version of A^-REPRESENTABILITY IILCV07I is also 
QMAEXp-complete: Given a density matrix p on two d-state particles, is it within e of a state p' such that 
there exists a translationally-invariant pure state [-(/') for particles arranged in a circle for which p' is the 
marginal state of two adjacent particles? p is a parameter of the problem, and N, given in binary, is the 
only input, as in our Hamiltonian problem. We can reduce to this version of A^-REPRESENTABILITY by 
starting with 1-DIM TIH on a circle. Then there is always a translationally-invaiiant pure ground state 
of the Hamiltonian H. By breaking the Hilbert space of two d-state particles up into small balls, we can 
get a finite set of density matrices p to try. For each one, if we can solve A^-REPRESENTABILITY, we can 
determine if p can be extended to a candidate ground state and if so we can determine the energy of\<j)), 
since it is just equal to Ntr{Hip). Trying all possible p, we can thus find the ground state energy of H, up 
to some e-dependent precision. 



3 Hardness of TILING 

The construction will make use of a binary counter Turing machine Mbc which starts with a blank semi- 
infinite tape. The head begins in a designated start state in the left-most position of the tape. Mbc will 
generate all binary strings in lexicographic order. More specifically, there is a function / : Z ^ {0, 1}* 
such that for some constant Nq and every N > Nq, if Mbc runs for N steps, then the string fsci^) will 
be written on the tape with the rest of the tape blank. Moreover there are constants ci and C2 such that if n 
is the length of the string fsciN) and N > Nq, then 2^^" < < 2=2" We will also assume that for any 
binary string x, we can compute A^ such that fsciN) = x'm time that is polynomial in the length of x. In 
some of the variations of the problem we consider in section [51 we will need to put additional restrictions 
on A^ (such as requiring A^ to be odd), and in those cases, we still require that we can find an A^ with the 
appropriate restrictions such that fsci^) = x. 

Using a standard padding argument, we can reduce any language in NEXP to NTIME(2'^i"). If L is in 
NTIME(2" ), the reduction consists of padding an input x so that its length is /ci |Pap95|. Thus, we 
will take an arbitrary non-deterministic Turing machine M which accepts a language L in time 2^^^" and 
reduce it to TILING. The tiling rules and boundary conditions will be specific to the Turing machine M but 
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Figure 1 : A possible tiling of the sides of a 5 x 5 grid. 

will be independent of any particular input. The reduction for Theorem 12.21 then will take an input string x 
and output integer N such that fBci^ — 3) = x. The tiling rules will have the property that a string x is in 
L if and only if an x grid can be tiled according to the tiling rules. 

Proof of Theorem 12.21 The boundary conditions for the N x N grid will be that the four corners of the 
grid must have a designated tile type (We will actually only need to use the upper left and bottom left 
comers.) First we will specify a set of boundary tiles and their constraints. In addition to ^ there are four 
other kinds of boundary tiles: 0, 0, We will call the rest of the tiles interior tiles. |^ will mark 
the left side of the grid, ^ the top of the grid, the bottom of the grid, and ^ the top side of the grid. 
(See figure [T]) 

Nothing can go to the left of a |^ tile which means that the only place a |^ tile could go is the left-most 
boundary, as desired. Similarly, nothing can above a tile, nothing can go below a tile, and nothing 
can go to the right of a ^ tile, which means @ tiles can only go in the top row, tiles can only go in 
the bottom row, and ^ tiles can only go in the right-most column. No interior tile can border a ^ tile in 
any direction. Furthermore a ^ cannot border on itself in any direction. This means that the only possible 
locations for a ^ are the four comers since those are the only places which can be surrounded by |^, 0, 
[ttttI , or ^ tiles. Since the boundary conditions state that ^ tiles must go in the corners, those are exactly 
the locations that will hold ^ tiles. The only tiles that can go above or below the ^ tiles are |^ and ^ 
tiles. The only tiles that can go either to the left or right of the ^ tiles are ^ and tiles. We will add the 
constraint that the only tiles that can go above or below tiles are |^ tiles or ^ tiles. Thus the entire west 
boundary, except for the comers, will be |^ tiles. Similarly, we add constraints that tiles must have ^ 
or tiles to their right and left and any tile must have a ^ or tile to its right and left. This means 
that the entire south border except for the corners is tiled with tiles and the entire north border except for 
the comers will be tiles. We do not need to put any further restrictions on the right boundary. 

The remainder of the grid will be tiled in two layers. The constraints on the two layers only interact at 
the bottom of the grid, so we describe each layer separately. The actual type for an interior tile is specified 
by a pair denoting its layer 1 type and layer 2 type. The bottom layer will be used to simulate the Turing 
machine Mbc- The top boundary of the grid will be used to ensure that Mbc begins with the proper initial 
conditions. Then the rules will enforce that each row of the tiling going downwards advances the Turing 
machine Mbc by one step. At the bottom of the grid, the output is copied onto layer 2. Layer 2 is then used 
to simulate a generic non-deterministic Turing machine on the input copied from Layer 1. The lower left 
corner is used to initialize the state of M and the constraints enforce that each row going upwards advances 
the Turing machine M by one step. Finally, the only states of M that are allowed to be below an tile are 
accepting states. Since each Turing machine only executes for A^ — 3 steps and the grid has space for A^ — 2 
tape symbols, the right end of the tape will never be reached. 

Although it is well known that tiling rules are Turing complete IIBer66l . we review the ideas here in 
order to specify the details in our construction. We will assume that the Turing machine M is encoded in a 
tiling that goes from bottom to top. This can easily be reversed for M bc which goes from top to bottom. 
The non-deterministic Turing machine M is specified by a triplet (S, Q, 5), with designated blank symbol 
# G S, start state qo ^ Q and accept state qa € Q- There are three varieties of tiles, designated by elements 
of S (variety 1), S x Q (variety 2) and S x Q x {R, L} (variety 3). Variety 1 represents the state of the 
tape away from the Turing machine head. Variety 2 represents the state of the tape and head when the head 
has moved on to a location but before it has acted. Variety 3 represents the state of the tape and head after 
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the head has acted, and the {R, L} symbol tells us which way the head moved. In the horizontal direction a 
tile corresponding to a G S can go next to any other interior tile to the left or right. Variety 2 types cannot 
go next to each other nor can variety 3 types. The only allowed pairings of variety 2 tiles and variety 3 tiles 
in the horizontal direction are of the form: [a, q] to the left of [6, g, L] or [a, q] to the right of [6, q, R]. Note 
that the state q of the head must be the same and the variety 2 tile must be placed in the direction designated 
by the i? or L in the variety 3 tile. 

In the vertical direction, for any a G S we can have [a] above [a]. Variety 2 tiles cannot go above or 
below each other and variety 3 tiles cannot go above or below each other. A variety 1 or a variety 2 tile can 
go above a variety 3 or a variety 1 tile as long as the alphabet symbols are the same. That is, [a, q] or [a] 
can go above either [a, q',L/R] or [a]. Finally, a variety 3 tile must go above a variety 2 tile and a variety 2 
tile must go below a variety 3 tile. Furthermore, these parings are allowed only if they encode a valid move 
of the Turing machine. That is, [b, q', L] can go above [a, q] only if (a, q) — > (6, q', L) is one of the valid 
non-deterministic moves of the machine. Similarly [5, q\ R] can go above [a, q] only if (a, q) {b, q' , R) 
is one of the valid non-deterministic moves of the machine. The table below gives an example of a section 
of tiles that encodes the move {a,q) —>■ {h,q',L) of the Turing machine: 



[c, q'] 


[b,q',L] 


[d] 


[c] 


[a,q] 


[d, q, L] 



The lower row shows the head in the square with the a. The [d, q, L] is from the previous TM move. The 
tile [b,q',L] enforces that the tiler is committing to executing the step {a,q) — > {b,q',L), although there 
may have been other non-deterministic choices. The [c,q'] tile to the left of the [b,q',L] shows the new 
location and state of the head. The [b, q', L] tile now just acts as a [b] tile for purposes of the tiling above. 

For our particular construction, we would like to start out the Turing machine Mbc with [go, #] in the 
leftmost location followed by [#] tiles. For layer 1, the only tiles that can go below a ^ tile are [qo, or 
[#] tiles. We forbid having a [#] tile to the right of a |^ tile and we do not allow a [qo , tile to the right 
of a [#] tile. We can assume without loss of generality that the Turing machine overwrites the leftmost # 
on the tape and never writes a # there again. We can also assume that the Turing machine never transitions 
back to the qo state. The rest of the layer 1 rules just enforce the rules for the Turing machine Mbc- 

Now in order to copy the output from Mbc to the input tape for M, we restrict the kinds of tiles that can 
go above tiles. A layer 2 tile that goes above a must be [a] or [a, qo] for some a G S. Furthermore, in 
the space above an tile, the alphabet characters for the layer 1 and layer 2 tiles must match. This copies 
the output of Mbc onto the input of V. Now we want to ensure that the starting configuration of V has only 
one head in the leftmost location. To accomplish this, we forbid a to go next to an [a] tile for a G S Mbc 
and forbid an [a] tile (for all a G S) to the left of a [b, qo] tile. Again, we can assume that M never transitions 
back to qo. A little care must be taken to overwrite the leftmost input tape character with something that is 
not in the alphabet of Mbc- This is because we have forbidden having an [a] tile to the right of a for 
any a G The information encoded in the left-most tape symbol can be retained by having a new a' 

symbol in S^/ for every a G '^Mbc- 

Finally, the only variety 2 tiles on layer 2 which we allow below a tile must be of the form [a, qA], 
where qA is the accepting state. Thus, there is a valid tiling if and only if the non-deterministic TM M 
accepts on input x in — 3 steps. 

4 The Quantum Case 

As in the 2-dimensional classical tiling problem, we make use of a binary counting Turing machine Mbc- 
Because we are working with quantum systems, we will require that M bc be reversible. Bernstein and 
Vazirani 1BV97] have shown that any deterministic Turing machine can be made reversible, meaning that 
given a configuration of the Turing machine, it has a unique predecessor in the computation. There may be 
some additional overhead but it is not significant. We can still assume that there is a function / : Z ^ {0, 1}* 
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such that for some constant A'o and every N > Nq, if Mbc runs for N steps, then the string fsci^) will 
be written on the tape with the rest of the tape blank. Moreover there are constants ci and C2 such that if n 
is the length of the string fsci^) ^rid A'' > Nq, then 2^^" < N < 2*^2" Wg also assume that for any 
binary string x, we can compute N such that fBci^) = a; in time that is polynomial in the length of x. 

We can reduce any language in QMA^xp to a language L that is accepted by a verifier who uses a 
witness of size 2*^1 " and whose computation lasts for 2'^^^ steps, where n is the length of the input. This is 
the same reduction used in the classical case, in which the input is padded to length \x\^/ci. We can use 
standard boosting techniques to assume that the probability of acceptance or rejection is at least 1 — e or at 
most e for e = l/poly(A^) [KSV02|. Suppose we are given an arbitrary verifier quantum Turing machine V 
which takes as input a classical/quantum pair {x, lip)) such that |V') has 2'^i" qubits and halts in 2'^i" steps. 
Based on V, we will produce a Hamiltonian term H which acts on a pair of finite-dimensional particles. We 
will also produce two polynomials p and q. The reduction for Theorem 12.51 will then take input string x 
and output an integer N such that fBci^ — 3) = x. The Hamiltonian will have the property that for any 
X, if there exists a \ip) that causes V to accept with probability at least 1 — e, then when H is applied to 
every neighboring pair in a chain of length A^, the resulting system has a unique ground state whose energy 
is at most p{N). If for every \ip), M accepts with probability at most e, then the ground state energy of the 
system is at least p(A^) + l/q{N). 

Quantum Turing machines were first introduced in IIDeu85ll and further developed in IIBV971 . The latter 
paper showed that we can make a number of simplifying assumptions about the form of a quantum Turing 
machine and not restrict its power in a complexity-theoretic sense. In particular, we can assume without loss 
of generality that the Turing machine V has a one-way infinite tape and that the head starts in designated 
start state qq at the left-most end of the tape. We will also assume that on input x, after 2^^^^^ steps, the 
Turing machine is in an accepting or rejecting state and the head is again at the left-most end of the tape. We 
will also assume that the witness will be stored in a parallel track with the left-most qubit in the left-most 
position of the tape. 

We now describe the set of states for the particles. A standard basis state for the whole system will be 
denoted by the state for each particle. States © and @ are special bracket states that occur at the ends of 
the chain. 

Definition 4.1 A standard basis state is bracketed ;/ the left-most particle is in state '<>, the right-most 
particle is in state > , and no other particle in the chain is in state ^ or Sbr is the space spanned by all 
bracketed states. 

We will restrict our attention for now to Shr and add a term later that gives an energy penalty to any state 
outside Shr- The rest of the states will be divided into six tracks, so the state of a particle is an ordered 6-tuple 
with each entry specifying the state for a particular track. The set of allowable states will not necessarily be 
the full cross product of the states for each track. 

Two of the tracks will implement a clock, with one track working as sort of a second hand and another 
track as a minute hand. The other four tracks will be used to implement two Turing machines which share a 
work tape. Track 3 holds the work tape. Track 4 holds the state and head location for the first Turing Machine 
(which is Mbc) and Track 5 holds the state and head location for the second Turing Machine (which is V). 
The sixth track will hold the quantum witness for V. Since there is limited interaction between the tracks, 
it will be simpler to describe the Hamiltonian as it acts on each track separately and then describe how they 
interact. The figure below gives a picture of the start state for the system. Each column represents the state 
of an individual particle. 
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Track 4: Tape head and state for TM Mbc 
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Track 5: Tape head and state for TM V 
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As is typical in hardness results for finding ground state energies, the Hamiltonian appUed to each pair 
will consist of a sum of terms of which there are two types. Type I terms will have the form \ab) {ab\ where a 
and b are possible states. This has the effect of adding an energy penalty to any state which has a particle in 
state a to the immediate left of a particle in state 5. We will say a configuration is legal if it does not violate 
any Type I constraints. Type IT terms will have the form: ^{\ab) {ab\ + \cd) {cd\ — \ab) {cd\ — \cd) {ab\). These 
terms enforce that for any eigenstate with zero energy, if there is a configuration A with two neighboring 
particles in states a and b, there must be a configuration B with equal amplitude that is the same as A except 
that a and b are replaced by c and d. Even though a Type 11 term is symmetric, we associate a direction 
with it by denoting it with ab cd. Type 11 terms are also referred to as transition rules. We will say 
that configuration A transitions into configuration B by rule ab cd if B can be obtained from A by 
replacing an occurrence of ab with an occurrence of cd. We say that the transition rule applies to A in the 
forward direction and applies to B in the backwards direction. We will choose the terms so that for any 
legal configuration, at most one transition rule applies to it in the forward direction and at most one rule 
applies in the backwards direction. Thus, a state satisfying all Type 1 and Type II constraints must consist of 
an equal superposition of legal configurations such that there is exactly one transition rule that carries each 
configuration to the next. The illegal pairs are chosen so that any state which satisfies the Type 1 and Type 11 
constraints corresponds to a process we would like to simulate or encode in the ground state. In our case, the 
process is the execution of two Turing Machines each for — 3 steps, where A'^ is the length of the chain. 

We will make use of the following simple lemma throughout the construction in Umiting the set of 
standard basis states in the support of the ground state. 

Lemma 4.2 For any regular expression over the set of particle states in which each state appears at most 
once, we can use illegal pairs to ensure that any legal standard basis state for the system is a substring of a 
string in the regular set. 

Proof: The alphabet for the regular expression is the set of particle states. Since each character appears once 
in the regular expression, the set of characters b which can follow a particular character a is well defined and 
does not depend on where the character appears in the string. Therefore, we can add an illegal pair ab if b is 
not one of the characters which can follow a. Any substring of the regular expression has no illegal pairs. ■ 

We give now a brief outhne of the construction and provide the full details in the subsections that follow. 
Illegal pairs are used to enforce that the state of Track 1 is always of the form @© (@ + ©)0 ©• (The + 
denotes the regular expression OR and not a quantum superposition.) There is one arrow symbol on Track 
1 that shuttles back and forth between the left end and the right end and operates as a second hand for our 
clock. We call one round trip of the arrow on Track 1 an iteration. Every iteration has 2(A'^ — 2) distinct 
states and 2{N — 2) transitions. Each iteration causes one change in the configuration on Track 2 which acts 
then as a minute hand for the clock. The Track 2 states are partitioned into two phases. The first phase is 
called the Counting Phase and consists of all iV — 2 of the states of the form 1) ©(Q) ©. The second 
phase is the Computation Phase and consists of all — 2 of the states of the form '^T*^ 2 *© states. 
The 1^(1;*®®*© states are ordered according to the number of particles in state ® and the v^®*®@*© 
states are ordered according to the number of particles in state @. The state immediately after ©® ®@ in 
the ordering is ® © ©@. The target ground state for the clock is the uniform superposition of all the clock 
states, entangled appropriately with states of the other 4 tracks. We need to have illegal pairs that cause all 
other states to have an energy cost. As is the case in other such proofs, it is not possible to disallow all states 
directly with illegal pairs. Instead, we need to show that some states are unfavorable because they evolve 
via forward or backwards transitions to high energy states. 

Each of the arrow states for Track 1 will come in three varieties: ^ and © will be used during the initial 
minute of the clock when it is in state ^ 'o"U^*'5' and will be used to check initial conditions on the other 
tracks. <t/ and ir^ will be used during the counting phase and <i/ and ir^ will be used during the computation 
phase. © and © will be used to trigger different actions on the other tapes. Every time the © sweeps from 
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the left end of the chain to the right end of the chain, it causes Mbc to execute one more step. Thus, Mbc is 
run for exactly N — 2 steps. The @ symbol is what causes the Turing machine V to execute a step. We then 
add a term that penalizes any state which is in the final clock state and does not have an accepting Turing 
machine state. Thus, only accepting computations will have low energy. 

Finally we use an additional term to enforce the boundary conditions. This is achieved by weighting 
the Hamiltonian terms for the illegal pairs and transition rules by a factor of three. Then an additional term 
is applied to each particle, which gives a benefit to any particle that is in the ? or state. Only the left- 
most and right-most particles can obtain this energy benefit without incurring the higher cost of having an 
endpoint state in the middle of the chain. 

4.1 The Clock Tracks 

There will be four types of states in Track 1: @, ©, ©, and O- (However, the @ and © states each come in 
multiple varieties — see below for the details.) We will have illegal pairs that will enforce that any standard 
basis state with no illegal pairs must be a substring of a string from the regular expression: 



Furthermore, any bracketed standard basis state must be exactly a string corresponding to this regular ex- 
pression. The transition rules then move the arrows in the direction in which they are pointing: 

1. @0 —>■ ©@5 ©© ©O: in the forward direction, arrows move in the direction they are pointing. 
In the reverse direction, they move in the opposite direction. 

2. @@ —>■ ©®, ®© — s- @@: arrows change direction when they hit an endpoint. 

Note that the symbol always interacts with the particle on its right in the forward direction and the 
particle on its left in the reverse direction. Similarly, the <- symbol interacts with the particle on its left in 
the forward direction and the particle on the right in the reverse direction. Therefore, we know that every 
configuration for Track 1 has exactly one transition rule that applies in the forward direction and one that 
appUes in the reverse direction. The arrow symbol on Track 1 shuttles back and forth between the left end 
and the right end as shown below and operates as a second hand for our clock. We call one round trip of the 
arrow on Track 1 an iteration. Every iteration has 2(A^ — 2) distinct states and 2(A^ — 2) transitions. Each 
iteration causes one change in the configuration on Track 2 which acts then as a minute hand for the clock. 



There are five states for Track 2: ®, ®, (2)and ©and ® . We will have illegal pairs that will enforce that 
any legal configuration on Track 2 must have the form @®*((Q)(fi)* + ®(2/*){>'. We will impose an ordering 
on the set of all possible such states and select transitions so that in one Track 1 iteration, the state on Track 
2 will advance from one configuration to the next in the ordering. Each state for Track 2 corresponds to a 
minute in our clock and the states are partitioned into two phases. The first phase is called the Counting 
Phase and consists of all the states of the form ©® (q)® ®- The second phase is the Computation Phase 
and consists of all of the states of the form '^® ®@ @. The '<"T' 'o"0^ > states are ordered according 
to the number of particles in state i and the '?*®*®@ (5^ states are ordered according to the number of 
particles in state The state immediately after (<)®*(q)c^ in the ordering ii, The ordering of 

the states for Track 2 is shown below for a six particle system. 
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©®®®®@ 
®®®®®@ 

The arrow particle on Track 1 will cause a transition from one Track 2 state to the next. Each of the 
arrow states for Track 1 will come in three varieties: t and V will be used during the initial minute of the 
clock when it is in state ©@® @ and will be used to check initial conditions on the other tracks. © and © 
will be used during the counting phase and @ and © will be used during the computation phase. We need 
to describe how these different Track 1 symbols trigger transitions on Track 2. The transitions for Track 1 
will remain as described in that the transition ^j, — > t; will always cause a transition from one of the 
@ symbols to one of the other Q) symbols, but we need to specify which @ symbols are used in order to 
ensure that forward and backwards transitions remain unique for each standard basis state. 

We will use ordered pairs to denote a combined Track 1 and Track 2 state for a particle as in [©,©]. 
The states for the left-most and right-most particles are not divided into tracks and are simply < or (>). As 
space permits we will denote the states for the different tracks vertically aligned. Examples for neighboring 
particle states are given below: 
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We will use the symbol • to denote a variable state. Thus [@, ©] is used to denote any state in which the 
Track 1 state is -> . 

Clock States for the Initialization: When Track 1 has a @ or © symbol, we want to ensure that Track 
2 is in state @®®*@. Therefore, we disallow [©, ©] and [©, ©] for any ® that is not ® or ®. We have 
the usual transitions that advance the Track 1 arrow: 
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These happen regardless of the values on the other tracks and do not change the values on the other Tracks. 
The state ©[©,(§.] does not have a transition in the reverse direction since this only occurs in the initial 
state. Finally, we have the transition 
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The presence of the 'o on Track 2 specifies that @ should transition to © in the reverse direction instead of 
©. The initial iteration of the second hand is demonstrated on a chain of length six in figure |2l 

Clock States for the Counting Phase: During the counting phase. Track 2 will always be in some 
state of the form < i o >, so we will forbid the states vi^] and [ir,, when • is 2 or i . The 
right-moving transitions will remain unchanged: -f — > ©@. These occur regardless of the contents of 
Track 2 and do not effect any change on Track 2. The change in direction at the right end will depend on the 
contents of Track 2: 
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The latter transition triggers the transition to the computation phase. In the left-moving direction, the arrow 
will sweep past pairs of ® particles and pairs of i particles: 
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Figure 2: The initial iteration of the second hand on a chain of length six 



When the anow on Track 1 sweeps left and meets the © particle on Track 2, it triggers an advance of the 
minute hand. This does not change the transition on Track 1 : 
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Note that the <r never coincides with the o , so we will disallow the state 
must turn at the left end: 



] . Finally, we have that © 





© 


® 


® 



© 



® 



We illustrate in figure [3] the first iteration of the counting phase and then the last iteration in the counting 
phase. The very last transition illustrates the transition to the computation phase. 

Clock States for the Computation Phase: During this time, Track 2 will always be in some state of the 
form ®® ®@ @, so we will forbid the states [©, ®] and [©, ®] when is ® or (a . Since should 
never be in the same configuration as ©, we will also disallow the pair [@, ®][0, ©]. The left-moving 
transitions will remain unchanged: '^ '^ . These occur regardless of the contents of Track 2 and do 

not effect any change on Track 2. The change in direction at the left end will happen as long as there is a ® 
left in Track 2: 



When the state on Track 2 becomes @®@ @ and the left-moving arrow on Track 1 reaches the left end of 
the chain, we want the clock to stop since this is the very last state. Thus, there is no forward transition out 
of ©[©, ®]. In the right-moving direction, the arrow will sweep past pairs of i particles and pairs of (2) 
particles: 
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When the arrow on Track 1 sweeps right and meets the ® particle on Track 2, it triggers an advance of the 
minute hand: 
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First iteration of the counting phase Last iteration of the counting phase 
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Figure 3: The first iteration and last iterations of the counting phase. The very last transition illustrates the 
transition to the computation phase. 
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Note that since the @ never coincides with the ®, we can make the state [@, ©] illegal. Finally, we have 
the turning at the right end: 
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We illustrate in figure |4] the first iteration in the computation phase and then the last iteration. The very last 
state shown is the final state for the clock. 

First iteration of the Computation Phase Last iteration of the Computation Phase 
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Figure 4: The first and last iterations of the computation phase. The very last state shown is the final state 
for the clock. 



Before describing the states and transitions for the other tracks, we will pause to consider the Hamil- 
tonian created so far. A single two-particle term will be the sum of the terms from the transition rules and 
the illegal pairs described so far. Suppose that a state of a particle is described only by its state on Tracks 
1 and 2. Let Hn be the system resulting from applying this term to every neighboring pair in a chain of 
particles. What is the dimension and structure of the ground state of Hj\[ restricted to S^r'^ 

Definition 4.3 We will say that a standard basis state is well-formed ;/ the Track 1 state corresponds to a 
string in the regular language vS^^i^^ (^ + ©)0 ^ and if the Track 2 state corresponds to a string in the 
regular language @®*(®®* -|- ®@*)@. 
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Note that if a standard basis state is well-formed it must be in Sbr- 



Lemma 4.4 Consider a well-formed standard basis state. There is at most one transition rule that applies 
to the state in the forward direction and at most one that applies to the state in the reverse direction. 
Furthermore, the set of well-formed states is closed under the transition rules. 

Proof: The rules are summarized in table |2l If only the Track 1 states are specified then the rule holds for 
any states on Track 2 and does not alter the Track 2 states. No pair appears twice on the left-hand side of a 
transition rule and no pair appears twice in the right-hand side of a transition rule. It can be verified that the 
application of a transition rule maintains the condition of being well-formed. 
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Table 2: The transition rules for Tracks 1 and 2. 



Consider the set of states that correspond to the valid clock states beginning with ®@0 ©on Track 
1 and ©®®*® on Track 2 and ending with @©0*@ on Track 1 and @®@*(>) on Track 2. There are 
exactly 4(A'^ — 2)^ such states in this sequence. Let be the uniform superposition of these states. 
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Lemma 4.5 Consider a well-formed standard basis state s that is not in the support of\(j)ci)- Then for some 
i < 2N, s will reach a state with an illegal pair after i applications of the transition rules in either the 
forward or reverse direction. 

Proof: We will argue that for each well-formed standard basis state, it is either in the support of \(pd), has 
an illegal pair or is within 2N transitions of a state which contains an illegal pair. In some cases, we have 
added penalties for a particular particle state (e.g., [9 , 2 ]). This can easily be handled with illegal pairs by 
making any pair illegal which contains that particular state. 

Every Track 1 configuration of the form @@ (@ + ©)0 ® occurs in \(j)ci)- Moreover, for every 
standard basis state in the support of \(j)ci), if the control state is @ or ©, then Track 2 is in state @(s)(Q)*@. 
Now consider a state in which Track 2 has a i , i or 2 particle and the control particle on Track 1 is @ or 
V . The control particle will transition (in either the forward or reverse direction) towards the X'-, (i) or 
particle on Track 2 and eventually they will coincide. This happens in fewer than moves and will result 
in an illegal state. 

Now consider the standard basis states whose control particle is @ or ©. Every possible combination of 
®©*(^ + ^ ) * > occurs in \(pci) with every possible combination of '^T* o * > except those where 
the control particle is *r and coincides with the o . The state [ *r , o ] is an illegal particle state. So we now 
need to take care of the case where the control particle on Track 1 is @ or © and Track 2 is a @®*®(D*@ 
state. Since it is illegal for a © or © to coincide with a i or 2 on Track 2, we can assume that the control 
particle on Track 1 is over a ^ particle on Track 2. In this case, it will transition in the forward direction in 
fewer than 2N steps (possibly turning at the left end of the chain) until it reaches the iX*- This will result in 
a particle in state [@, ®] which is illegal. 

Finally, consider the case where the control particle is 9^ or ©. Every possible combination of the 
expression + ©)0*@ occurs in \(/)ci) with every possible combination of @®*(X'(2i* > except 

those where the control particle is and coincides with the X- The state , i ] is illegal. So we now 
need to take care of the case where the control particle on Track I is © or © and Track 2 is a 
State. Since it is illegal for a © or © to coincide with a (Q) or ® on Track 2, we can assume that the control 
particle on Track 1 is over a T particle on Track 2. In this case, it will transition in the forward direction in 
fewer than 2N steps (possibly turning at the left end of the chain) until it is just before the 6 on Track 2. 
This will result in a pair [©, ®] [©, ©] which is illegal. 

■ 

Lemma 4.6 \4>d) the unique ground state of H]^\s^^. All other eigenstates have an energy that is at least 

n{i/N^). 

Proof: The proof of this lemma follows from the standard techniques used for showing QMA-completeness 
of 1-dimensional Hamiltonians [AGIK07I . so we only give a brief sketch here. The idea is that any standard 
basis state inside Sbr which is not well-formed will have energy at least one from the illegal pairs, so we 
can restrict our attention to well-formed states. Now we create a graph over well-formed standard basis 
states. There is an edge from state a to state 5 if 6 can be reached from a by the application of one transition 
rule in the forward direction. Lemma l4!4l implies that this graph is composed of disjoint directed paths. 
Call Htrans the Hamiltonian resulting from the sum of all the terms from transition rules and Hi^gai the 
Hamiltonian from illegal pairs. The subspace spanned by the states in a single path is closed under Htrans 
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and Hiegai- Furthermore, the matrix for H trans restricted to the states in a path looks like 
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where the dimension of the matrix is the length of the path. The matrix for Hi^gai is diagonal in the standard 
basis. If the path contains no states with illegal pairs, then ///ega/ is zero. The unique ground state has zero 
energy and is the uniform superposition of all standard basis states in the path. Moreover, the next highest 
energy state has energy at least 1/L, where L is the length of the path. In our case. Lemma 1431 indicates 
there is exactly one path with no illegal states which corresponds to \4>ci)- The length of this path is 0{N'^), 
so the next highest eigenvalue in this space is at least Q{1/N'^). 

Now consider a path with some illegal states and suppose that the ratio of illegal states in the path to the 
length of the path is 1/s. Again, it is known IIKSV02II that any state in the subspace spanned by the states 
in this path will have energy be at least 1/s^. By Lemma 1431 we know that for any path which does not 
con^espond to \(pci) the ratio of illegal states to the total number of states at least 1/2N, which means that 
any state in this subspace will have energy at least ^l{l/N^). ■ 

The remaining subsections describe how the arrow on Track 1 interacts with the other four tracks. The 
only role of Track 2 is to record the time. It causes the control state on Track 1 to transition from the counting 
phase to the computation phase and finally to stop iterating at the end of the computation phase. 

4.2 Initialization Phase 

Track 3 will be used as the tape for both Turing machines. Therefore the set of track 3 states for a particle 
consists of the union of the two alphabets for the two Turing machines. We want to ensure that this track is 
initialized to all blank symbols, so we disallow any state in which Track 1 is in and Track 3 has anything 
butC- 

Track 4 will store the location and state for Mbc- Therefore, we will add illegal pairs to ensure that the 
Track 4 state for the system always has the form @@ ©O C> , where • is a state for Mbc- We assume 
that Mbc has a designated start state S , so we want to ensure that Track 4 starts out in the configuration 
©©O*®- To do this, we make any state illegal in which the Track 1 state is @ and the Track 4 state is 
anything except (§) or O- Furthermore, for any pair in which the left particle is in state @, and the Track 
1 state is @, the Track 4 state must be S). For any pair in which the left particle is not in state @, and the 
Track 1 state is ? , the Track 4 state must be ^ . The initial conditions for Track 5 are similar, except that 
the starting state for V is used. Since Track 6 holds the quantum witness for V, it can be any state on n 
qubits. A state of a particle is now specified by a 6-tuple. 

4.3 Counting Phase 

Every time the © sweeps from the left end of the chain to the right end of the chain, it causes Mbc to 
execute one more step. Thus, Mbc is run for exactly N — 2 steps. Since the classical, reversible Turing 
machine Mbc is a special case of a quantum Turing machine, we describe the details for the more general 
case. The @ symbol is what causes the Turing machine V to execute a step. 
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4.4 Computation Phase 



We will examine a particular quantum rule and explain the desired behavior of our machine. Consider a pair 
{q, a), where q ^ Q and a G S. a will encompass the state on the work tape as well as the state on the 
witness tape. Since the @ particle triggers the execution of a step, we will consider particles as triplet states 
of the form ['■? , g, a], ['"'' , q, a] and ["^ , g, a], for a € S and q e Q U {O} U { ; • }. q specifies the Track 5 
state, and a tells us the state of Track 3 and Track 6. It will be convenient to refer to g as a generic Track 
5 state of a particle which could be a state from Q as well as O or © . We will also be interested in the 
computation state of a particle which just consists of a pair [q, a], where q ^ Q U {O} U {©} and a € H. 

If the TM is in state q and in a location with an a on the tape, the QTM defines a superposition of 
possible next moves. Let 6{q, a,p, b, D) denote the amplitude that the next state will be a configuration in 
which the state is p, the tape symbol is overwritten by b and the head moves in direction D. We will use a 
fact established in IIBV97II that the states of a quantum Turing machine Q can be partitioned into two sets 
Ql and Qr such that states in Ql can be reached only by moves in which the head moves left and Qr can 
be reached only by moves in which the head moves right. We will use q^ to designate a generic element of 
Ql and q^ to designate a generic element of Qr. 

We will need to execute the moves in which the head moves left separately from the moves in which the 
head moves right. In order to do this, we introduce a new state q'^ for every state qji G Qr. We will call 
this set Q'j^. It will be forbidden to have the Turing machine head in a state from Q'^ without also having 
the @ state on Track 1 in the same location. It will also be forbidden to have the (?) in the same location as 
g if g € Ql- We will be interested in two particular subspaces defined on the computation state for a pair of 
neighboring particles [q, a][p, b], where a, 6 € S and q,p ^ Q U {O} U {©}. Sa is the space spanned by 
all such legal basis states such that q Ql, p Q'r- Note that if the Track 1 state for a pair of particles is 
|@0), then the only possible computation states for the pair are in Sa- We also define Sb to be the space 
spanned by all legal basis states of the form [q, a][p, b] such that p ^ Ql, Q ^ Q'r- Similarly, if the Track 1 
state for a pair of particles is | ? ), then the only possible legal computation states for the pair are in Sb- 

We will describe a transformation T that maps Sa to Sb- The extended map T (g) | )(:^0| as it 
is applied to each pair of particles in turn will implement a step of the Turing machine. We will call this 
extended map a transition rule since it carries one state to another in the same way that the transition rules 
for the clock state did. The difference now is that the states will in general be quantum states. 

The transformation on Sa works in two parts. At the moment that the Track 1 arrow moves from the 
position to the left of the head, we execute the move for that location. For every {qL, b, L), with amplitude 
6{q,a,qL,b,L), we write a 6 in the old head location and move the head left into state qL- For every 
{qn, b, R), with amplitude 6{q, a, q^, b, R), we write a 6 in the old head location, transition to state q'p, and 
leave the head in the same location. We need to defer the action of moving the head right until we have 
access to the new location. In the next step of the clock, when the ^ is aligned with the state q'j^, we move 
it right and convert it to q^. 

That is, we want a sequence of two transitions. 



© 



^ S{q,a,qL,b,L) 



+ X] ^i9^a,qR,b,R) 
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(1) 



After this step, the configuration is a superposition of states in which the step has been performed on the 
configurations with the head in the same location as the -f state, except that moving the head to the right has 
been deferred. If the TM state is primed and is aligned with the location of the , that triggers the execution 
of the right-moving step. 
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Otherwise, the @ state just sweeps to the right, leaving the other tracks unchanged: 
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(3) 



Note that transformations ([T|), Q and ^ define T on every state in Sa- Furthermore, the image of T\s^ is 
in Sb- a critical fact is that after T is applied to a pair, the computation state of that pair is in Sb, which 
implies that the next pair over is now in subspace Sa and T can be applied to this new pair. We need to 
now show that T has the necessary properties to be expressed as a Type II Hamiltonian term with the \cd) 
final state a superposition. That is, we want to show that T is a partial isometry (meaning it preserves inner 
products) when restricted to Sa- Since Sa and Sb have the same dimension, this means that is well- 
defined and a partial isometry on 5^. We can extend T to be a unitary map on the full Hilbert space. Then 
we select the Hamiltonian term to be: 

Isa » I i'0)(@0| + Isb ® -T^ |TO)(@©| - Tt |©©)(@@|. 

Expressions ([T]), Q and ^ each define T on a different orthogonal subspace of Sa- Sai is the space of 
the form [@, C!^] [<?, a], Sa2 is the space of states of the form [g^, ©] [©, ®] and 5^3 is the space spanned 
by states of the form [qr, ©][©, ©], [©, ©][©, ©] and [©, ©][©, ©]. We will argue that T is a partial 
isometry when restricted to each of Sai, Sa2 and Sa3- Furthermore the images of Tg^i , ?5a2 ^^'^ '^Sa3 
mutually orthogonal. This implies that T is a partial isometry on Sa- 

T is clearly a partial isometry on each of 5^2 and 5^3- We only need to examine the transformation 
given in ([T]). We know that the application of a move of the QTM is unitary. In particular, it is unitary 
when restricted to the space of all configurations with the head in a particular location. This means that the 
transformation ([T]) (T restricted to states of the form [ ,•][(/, a]) is a partial isometry. 

Finally, we must compare Sa in different locations. In particular, when the TM head moves right from 
location i, it can end up in the same location as when it moves left from location i+2. In our implementation, 
when each step is implemented, these steps are orthogonal because the Track 1 state differs. They are also 
orthogonal in the original TM rule, since after the head moves right, it will be in a state from Qr, whereas 
after it moves left, it will be in a state from Ql, which is orthogonal to Qr- Thus, our implementation is 
globally performing the correct TM rule. 

We have argued that the extended transformations (or transition rules) which implement the Turing 
Machine behave as we would like them to as long as the states to which they are applied obey certain 
conditions. The following definition formalizes those conditions. 

Definition 4.7 We say that a state is invalid if one of the following conditions holds: 

Furthermore, for some particle other 



1. The control particle on Track 1 is on site i and is in state t 
than i, its Track 4 state is q'^for some q'j^ € Q'j^. 

2. The control particle on Track 1 is on site i and is in state 
particle i is for some ql £ Ql. 

3. The control particle on Track 1 is on site i and is in state i 
than i, its Track 5 state is q'^for some q'j^ € Q'j^. 

4. The control particle on Track 1 is on site i and is in state 
particle i is ql for some qL ^ Ql- 



■ . Furthermore, the Track 4 state for 
Furthermore, for some particle other 
■ . Furthermore, the Track 5 state for 



We say that a standard basis state is valid if it is not invalid. Any invalid state has an energy penalty 
of at least one by the illegal pairs added in this subsection. The subspace of all valid states is closed under 
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the action of the transition rules. As before, we call a state well-formed if its clock state (Track 1 and 2) 
is well-formed. For any valid, well-formed state, the transition rules apply non-trivially to only one pair of 
particles. Furthermore, the transition rules are reversible and norm-preserving over the space spanned by 
valid, well-formed states. 

4.5 Combining the Tracks 

We have described a set of transition rules for Tracks 1 and 2 which advance a clock through T = 4(A^ — 2)^ 
states. We have also described a set of transition rules which implement two Turing Machines on tracks 3 
through 6. In each set of transition rules, the control particle of the clock on Track 1 advances one space. 
Thus, when they are combined, each transition rule advances the clock and implements part of a step on one 
of the Turing Machines. Consider any state such that the states for Tracks 1 and 2 are clock states and the 
states for Tracks 3 through 6 are a quantum state within the subspace of valid states. If the clock state is not 
the final state, then the transition rules apply to exactly one pair of particles and take the state to a unique 
new state with the same norm in which the clock has advanced by one tick. Similarly, if the clock state is not 
the initial state, then the transition rules applied in the reverse direction apply to exactly one pair of particles 
and take the state to a new state with the same norm in which the clock has gone back by one tick. 

We argued in Subsection 14.11 that a zero eigenstate of the Hamiltonian must be a superposition of the 
sequence of clock states. Now that we are considering the states of the other tracks as well, this is a higher- 
dimensional space. In order to define a basis for the zero eigenspace of the Hamiltonian defined so far, we 
specify a standard basis state for Tracks 3 through 6 for the initial clock state. It will be convenient to separate 
the states of Tracks 3 through 5 from the state of Track 6. So let X be a standard basis state for Tracks 3 
through 5 and 1" be a standai^d basis state for Track 6. Let \t) be the state of the clock after t transitions. 
Define denote the state of Tracks 3 through 6 after t applications of the transition rules assuming 

that Tracks 3 through 5 start in state X and Track 6 starts in state Y . Define |(/>x,y) = Ylt=o \^)\4't,x,Y)- 
The \4>x,y) satisfy all the constraints due to transition rules and illegal states for Tracks 1 and 2. We will 
use / to denote the desired initial configuration for Tracks 3 through 5. That is, Q on Track 3 and S;U 
and qo * on Tracks 4 and 5. The energy for any |0x,y) where X ^ I will be at least Q{1/N'^). This is 
because if X ^ I, at least one state in the initialization phase will have an energy penalty as the Q) control 
state sweeps to the right. So the subspace spanned by {|0/,y)} for all Y form a basis of the zero eigenspace 
for the Hamiltonian terms defined so far, restricted to S^r, the subspace of all bracketed states. The next 
highest eigenstate in Shr has energy 0,(1 /N^). 

4.6 Boundary Conditions 

First we address the case of a finite chain and add a term that will enforce that the ground state is in 
Sbr- In order to penalize states that are not bracketed, we weight H by a. factor of 3 and add in the term 
I - |@)(@| — to every particle. This has the effect of adding energy — 2 to every state in 5^^. 

The ground space is still spanned by the [0/,y) and the spectral gap remains Q{l/N^). Any standard basis 
state outside Shr will have a cost of at least — 1: If there are a particles in state >, or ^ in the middle 
of the chain, they will save a from the / — |@)(@| — |©)(©| term. However, there will be at least \a/2] 
illegal pairs because the @ particles in the middle will each form an illegal pair with the particle to their 
left and the > particles in the middle will each form an illegal pair with the particles to their right. We can 
lower bound the additional cost by counting the type ( < or > ) which is more plentiful, yielding at least 
[a/2] distinct illegal pairs. Thus, there will be a net additional energy cost of at least a/2 when a is even, 
or (a + 3)/2 when a is odd. 

If we instead have periodic boundary conditions, the problem is still QMAEXp-complete. We start by 
removing the penalty for the pair >, < . The set of legal and well-formed states is exactly the same as it was 
for the finite chain except that we can now have more than one segment around the cycle. For example, we 
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could have the following state wrapped around a cycle: 

© • • • @© ■ ■ ■ ©p ■ ■ • @ . 

Segment i Segment 2 Segment 3 

Note that the transition rules do not change the locations or numbers of @ or @ sites, so we can restrict our 
attention to subspaces in which the segments are fixed. If a standard basis state is well-formed then every 
occurrence of @ has a @ to its immediate right and every occurrence of © has a @ it its immediate left. 
Thus, we can assume that a standard basis state in the support of a ground state can be divided into valid 
segments. Of course, it is possible that there are no segments in which case the state could simply be a single 
particle state repeated around the entire cycle. We know that the states within a segment must be ground 
states for a chain of that length. Otherwise, the energy of the state is at least where I is the length 

of the segment. We need to add additional terms and states which give an energy penalty to states with no 
segments or more than one segment. 

It is most straightforward to achieve this by using odd N and adding an additional Track 0. There will 
be no transition rules that apply to Track 0, so the state of the Track remains fixed even as transition rules 
apply to the state. We will add some extra terms which are diagonal in the standard basis for Track and 
introduce energy penalties for certain pairs of particles. Fix a standard basis state for Track and a set of 
segments and consider the space spanned by these states. The Hamiltonian is closed over this space, so 
we just need to examine the eigenstates restricted to each such subspace. Track will have three possible 
states ®, ®, and Even the @ or @ particles have a Track state, so the state of a particle is either in 
{© , ® , CD} X {©, ©} or it is a 7-tuple denoting states for Tracks through 6. 

We add the terms [ A A ) ( A A ] and | b" b ) ( b b | , which gives an energy penalty of 1 for pairs A A or 
(B I3 on Track 0. Since is odd, the state will have energy at least one unless there is at least one particle 
whose Track state is ^L. Now we add an energy penalty of 1 for any particle whose Track state is J) 
which is not [©,©]. We also add an energy penalty for [®, ©] or [©,©]. The only way for a state to 
have energy less than one is for there to be at least one segment and to have a © on Track exactly at the 
locations where there is a '< on the other tracks. Finally, we add an energy penalty of 1/2 for any particle 
whose Track state is Uy. Thus, the only way for a state to have energy less than one is for there to be 
exactly one segment and to have a © on Track exactly at the single location where there is a © on the 
other tracks. 

Consider the set of all basis states with exactly one segment and whose Track state has a T co-located 
with the site and alternating vA/'s and vg^'s elsewhere on Track 0. Any eigenstate outside this space 
has energy at least one. Note that there are 2N different low-energy ways to arrange the Track states, 
since there are possible locations for the © site and the string of alternating 'A^'s and ^'s can begin 
with or 'B" . Each such choice gives rise to a subspace that is closed over H and otherwise has identical 
behavior with respect to the rest of the terms, so we will make an arbitrary choice and examine the resulting 
subspace. Since there is one segment, the state \4>i,y) from the previous subsection is well defined. This 
state is a ground state for H and has energy 1/2. All other eigenstates within the subspace have energy at 
least l/2 + n{l/N^). 

4.7 Accepting States 

Finally, we add a term for those 10/,^) that do not lead to an accepting computation for V. In the last 
configuration of the clock, the two right-most particles are in state [ < ] [ v , 1 ] . We assume that an accepting 
computation will end with the head in the left-most place in state qa- We add a penalty for a pair whose 
Track 1 and 2 state is [< ][<?, 1 ] such that the right particle in the pair does not have a Track 5 state of qa- 
Thus, an accepting computation will have an energy penalty of at most e/N"^ for this term where e can be 
made arbitrarily small. If no witness leads to an accepting computation, all states will have energy at least 
(1 — e) /N^ from this term. 
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Figure 5: Converting a boundary condition on two opposite comers to one on two adjacent corners. 

5 Variants of the Classical Tiling Problem 

We now turn to studying variants of TILING. There are variety of different boundary conditions we could 
consider. We can also consider changing the absolute prohibitions on certain adjacent tiles to a soft condition 
by assigning different weights to the different adjacent pairs of tiles. We can add additional reflection or 
rotation symmetry to the tiling rules. We do not have definitive results on the complexity of all combinations 
of these variations, but we have been able to determine the complexity of most combinations, and in this 
section we present our analysis. Also, we prove that in one dimension, TILING and all the above variants 
are easy. 



5.1 Boundary Conditions 

Our choice of fixing the tiles in all four comers of the square is unusual. See, for instance, Papadim- 



itriou |Pap95|, who fixes just the tile in a single corner. If we had instead chosen that convention, the 



TILING problem would become easy, in an annoying non-constructive sense: 

Theorem 5.1 Define a variant of TILING where a designated tile ti is placed in the upper right comer of 
the grid, and the other comers are unconstrained. Then there exists Nq E Z"*" U {oo} such that if N < Nq, 
there exists a valid tiling, and if N > Nq, then there does not exist a tiling. However, Nq is uncomputable 
as a function of (T, V). 

This theorem follows immediately from the observation that with this boundary condition, a valid tiling 
for an X grid can be cropped by removing the leftmost column and bottommost row to give a valid 
tiling for the (A^ — 1) x (A^ — 1) grid. We know Nq must be uncomputable because the question of whether 
there is an infinite tiling is uncomputable. Still, if we fix (T, H, V), we know there exists a straightforward 
algorithm to solve this variant of tiling: simply determine if A^ < Nq. We just do not know Nq, so we do 
not know precisely what algorithm to use. 

On the other hand, if we fix boundary conditions in two comers, that is already enough for hardness. Our 
proof (section |3]l uses pre-determined tiles in two adjacent corners, but two opposite corners would suffice 
as well. Suppose they are the top right and bottom left corners, and ^ is the tile type fixed for the upper 
right and lower left comers. Then we can introduce three additional tile types and 0. We allow 

only to be to the right of and only or ^ to be to the right of 0. We also allow only ^ to be 
below ^ and only ^ or ^ to be below There is one more rule that says that ^ cannot be above 
m or to the left of 0. There are other pairs, such as I UtttH , which we can either forbid or allow; they 
will not occur regardless. Naturally, the other tiles used in the construction will have allowed and forbidden 
configurations between themselves as before. 

Thus, the bottom row will be all or ^ (except for the left end), and the right column will be all ^ 
or 9 (except for the top). Then the bottom right corner can only be fixing the tile there as well. Indeed, 
PI can only be in the bottom right comer, so the bottom row is all except for the corners, and the right 
column is all ^ except for the corners. See figure [5l 

Another interesting case is when we have periodic boundary conditions. That is, we consider the top row 
to be adjacent to the bottom row, and the leftmost column is adjacent to the rightmost column. Essentially, 
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we are tiling a torus. This case is particularly interesting because it is truly translationally invariant, unlike 
our usual formulation, where the boundaries break the translational invariance. 

Theorem 5.2 Define PERIODIC TILING as a variant of TILING with periodic boundary conditions on an 
N X N grid. PERIODIC TILING is NEXP -complete under an expected poly-time randomized reduction or 
a deterministic polyspace reduction. 

Since we are dealing with the class NEXP, even an exponential-time reduction is meaningful. For 
instance, if there is an algorithm to solve PERIODIC TILING which takes time poly(A^), then, as a conse- 
quence of theoremO EXP = NEXP. 

Proof: To prove this, we will show that, for appropriate N, we can introduce an effective horizontal and 
vertical border at some point inside the square. The actual location of the borders will not be specified at all, 
but we will choose conditions so that there is exactly one horizontal border and one vertical border. Those 
borders will then act like the usual edges of the grid for the standard 4-corners boundary condition. 

We will specialize to N which is an odd prime. We will add two additional layers of tiles over those for 
the usual TILING NEXP-completeness result (section O. The new layer 1 has 7 different possible types of 
tile: B' IH' ffl' □' H' H' and Q. The new layer 2 has 10 types of tile: @, 0, U), [|], Q, Q, 

QandQ 

The and ^ tiles will create the vertical and horizontal borders. | and Q are used to make 

sure that there is at least one of each kind of border, and the Q and Q tiles will create a diagonal line within 
the rectangle defined by the borders. Layer 2 is used to mark the directions next to the border and make sure 
the diagonal line goes from the upper left comer of the rectangle to the bottom right comer, which ensures 
that the rectangle is a square. When is prime, this means there can only be one horizontal and one vertical 
border. 

Layer 1: Q can only be adjacent horizontally to Q or and ^ can only be adjacent vertically to 
m or ^ can only be adjacent vertically to and can only be adjacent horizontally to These 
rules imply that if we have a ^ anywhere in a column, that column must only contain ^ and ^ tiles, 
and if we have a Q anywhere in a row, that row can only contain Q and ^ tiles. Thus, layer 1 contains 
some number of vertical and horizontal lines composed of ^ and Q tiles intersecting at ^ tiles. No two 
horizontal or vertical lines can be adjacent. 

The next set of rules for layer 1 enforce that the space between the lines must be filled by a checkerboard 
of Q or altemating with | or Q tiles. We ensure this by forbidding Q and Q from being adjacent 
to themselves or each other in any direction, and forbidding | and Q from being adjacent to themselves 
or each other in any direction. Since N is odd, this ensures that it is not possible to tile the entire torus with 
the checkerboard pattern and there must be at least one horizontal and one vertical line. 

Furthermore, if we ever have a Q tile anywhere, we want to be forced to have a diagonal line of Q 
tiles continuing to the upper left and lower right, with a diagonal line of Q tiles next to it (above and to the 
right), with both lines ending at a vertical or horizontal line. We enforce this by requiring that a Q tile must 
have Q both below it and to its left. Above a Q tile we can have only Q or and to the right of a Q 
tile, we must have either ^ or 

Layer 2: We wish ^ to mark the top side of the rectangles delineated by the layer 1 Q and ^ tiles. 
m will mark the right side, ^ the left side, and the bottom side, with and ^ marking the 

upper left, upper right, lower left, and lower right comers, respectively. Q will be inside the rectangle and 
P will go over the layer 1 Q and ^ tiles. 

can have ^ adjacent to it in either horizontal direction; it could also have ^ to the left and ^ to 
the right. A tile can have ^ adjacent in either horizontal direction, and could have Q to its left or Q 
to its right. ^ can have ^ adjacent to it in either vertical direction, or could have ^ above or ^ below. 
in can have |^ adjacent to it in either vertical direction, or could have ^ above or ^ below. ^ must 
have a to its right and a |^ below it. ^ must have a to its left and a ^ below it. ^ must have a 
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Figure 6: A possible arrangement of layers 1 (a.) and 2 (b.) for PERIODIC TILING. 



nn above it and a to the right of it. ^ must have a to its left and a ^ above it. Thus, if we have a 

anywhere, there must be a horizontal line of that either goes all the way around the grid (due to the 
periodic boundary conditions) or ends at ^ on the left and ^ on the right. If it ends, then below the Q 
tile is a vertical line of |^ tiles which end at a ^ tile, and below the ^ tile is a vertical line of ^ tiles 
which end at a ^ tile. The ^ and ^ tiles must be in the same row, and between them is a line of 
tiles. Following this line of logic for the other tiles, we find that layer 2 must be a collection of rectangles, 
horizontal stripes, and vertical stripes. 

Now we add rules that enforce that we have | on the outside of the rectangles and Q on the inside. 
That is, if any of the border tiles has a solid color on one edge, it must be adjacent on that side to a matching 
solid color tile. ^ must have Q above it and Q below it, and it is the other way around for below 
and Q above). |^ must have Q to the left and Q to the right, and ^ must have Q to the left and Q 
to the right. Q must have Q above and to the left of it. ^ must have Q above and to the right of it. Q 
must have below it and to its left. ^ must have Q below it and to the right of it. Q can be adjacent 
to anything but except as noted above. Q can be adjacent to anything but except as noted above. 
These rules imply that inside the rectangle delineated by the 0, 0, and |^ tiles are only Q tiles, and 
immediately outside it are Q tiles. Each rectangle is outlined by 0, 0, and tiles on the sides, and 
B- CB' '^nd g on the corners, and is full of Q inside. A vertical stripe has on its left and ^ on its 
right, and a horizontal stripe has above it and below it. Both vertical and horizontal stripes have only 

1 I inside them. The rectangles and stripes cannot be adjacent and are separated by Q tiles. 

Interactions between layers 1 and 2: A layer 2 tile must be in the same spot as a layer 1 
or and only a Q tile can be in those spots on layer 2. The ^ ^^id ^ layer 2 tiles must be in the same 
spot as a layer 1 Q tile, but a layer 1 Q tile could also be in the same spot as a layer 2 Q tile. A layer 1 
Q tile cannot be in the same spot as any other layer 2 tile. 

Once we look at the compatibility conditions between the two layers, we see that in fact layer 2 can only 
contain rectangles which must be lined up with the space between the horizontal and vertical lines on layer 
1. On layer 1, each rectangle delineated by the horizontal and vertical lines must have a Q in the upper left 
corner, which starts a diagonal line of Q tiles extending towards to the bottom right. Since it must end at 
the border of the rectangle, but a Q tile cannot be in the same spot as a layer 2 ^ or tile, the only place 
the diagonal line can end is at the lower right corner. Thus, the rectangle must actually be a square. (See 
figure [6] for an example of an allowed layer 1 tiling.) 

However, the only way for all the rectangles formed by the horizontal and vertical layer 1 lines to be 
squares is if the spacing between them is equal. They then form M x M squares arrayed in a x grid for 
a total of k'^ squares. It follows that N = k{M + 1). But when N is prime, then k must be 1. {M cannot 
be since the horizontal and vertical lines cannot be adjacent.) Thus, the only allowed tiling is to produce a 
single {N — 1) x {N — 1) square on layer 1. We can then consider some data layers of tiles which implement 
the rules from section |3] The layer 2 Q, and ^ tiles mark out the corner of the square, so we can 

put a condition on the data layers that enforce the corner boundary conditions on those locations. 

We do need to be careful of one aspect, however, since we are now restricted to a size of square which 
is 1 less than a prime number. For any input x, we need to choose a prime N . N should be not much bigger 
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than X (log N = poly(log x)), and it must be possible for the universal TM implemented by the data layers 
to deduce x in a reasonable time. It is known that there exists a constant 6 < 1 such that 



lim[^(y + y^)-7r(y)] = -^, (4) 
j/^oo log y 



where 7r(y) is the number of primes less than or equal to y | Ing37[ . (For our purposes, it is sufficient to 



take 9 = 2/3, but smaller 6 is possible.) That is, for sufficiently large y, the number of primes in the 
interval [y,y + y^] is approximately 1/logy times the size of the intervalQ Thus, given x, we will let 
no{x) = 2[logj;] (that is, basically twice the number of bits in x) and let No{x) = j;2"o(^) (that is, the 
binary expansion of No{x) is that of x followed by {no{x) — 1) Os). Then 

Nq(^x)^ = x^2^""^'*^ < 2^(1"^°^ ^"l+^o^^')) < 2""^"^^ (5) 

In particular, it follows that Nq{x) + A^o(a^)^ ^ Nq{x + 1). If we can find a prime in the interval 
I{x) = [N(){x), Nq{x) + NQ{xy), the TM can thus easily deduce x by looking at the most significant bits. 
There is at least one prime in I{x), so we can certainly find one with an exhaustive search, which can be 
done with polynomial space. Furthermore, by if we choose a random N G I{x), for sufficiently large x, 
there is a probability about l/log(A'^o(3;)) that N is prime, which we can verify in poly (log x) time. Thus, 
we get a randomized reduction to PERIODIC TILING which runs in expected time poly (log x). 



5.2 Weighted Constraints 

Now we consider the case where the constraints count different amounts. 



Definition 5.3 WEIGHTED TILING 

Problem Parameters: A set of tiles T = {ti, . . . , tm}- A set of horizontal weights wr : T x T — > Z, such 
that ifti is placed to the left oftj, there is a contribution ofwH{ti,tj) to the total cost of the tiling. A set of 
vertical weights wy : T x T ^ 7^, such that ifti placed below tj, there is a contribution ofwv{ti, tj) to 
the total cost of the tiling. An integer c. 
Problem Input: Integer N, specified in binary. 

Output: Determine whether there is a tiling ofanN x grid such that the total cost is at most c. 

Notice that in this case we do not make any restrictions on the boundary conditions. We can define an 
additional variant of WEIGHTED TILING where we instead have periodic boundary conditions (the top 
and bottom rows are considered adjacent, as are the left-most and right-most columns). 

The definition of WEIGHTED TILING is not completely analogous to the quantum 1-DIM TIH prob- 
lem. We require the maximum cost to be a constant independent of A^. More generally, we could allow c 
to be a polynomial in A'^. This does not change the complexity of any case for which we have been able to 
prove a definite result, with the exception of WEIGHTED TILING with reflection symmetry and periodic 
boundary conditions, where we have shown the problem is in P when c is constant, and have no result when 
c = VL{N). 

Theorem 5.4 WEIGHTED TILING is H^YS' -complete, as is the periodic boundary conditions version of 
WEIGHTED TILING. 

Proof: WEIGHTED TILING with open boundary conditions: It is sufficient to take all the weights to 
be 0, +2, or —1, and c = —4. We will reduce from TILING; the main point is to show we can fix the 
boundary conditions on the comers. To do so, we take the TILING instance and add a new layer of tiles 

' This is the result one might expect from the Prime Number Theorem, but that theorem is not strong enough, as it is compatible 
with having a large interval with a low density of primes. 



26 



1 1 


II 


II 


1 1 


1 1 


II 


II 


1 1 


1 1 


II 


II 


1 1 




II 


II 


II 



Figure 7: The preferred anangement of the new layer in WEIGHTED TILING with open boundary condi- 
tions. 

consisting of five special types of tile, and ^ will have a weight +2 if it is below any 

type of tile at all, and a weight +2 if it is to the left of any type of tile at all. However, it has a weight —1 if it 
is to the right of any type of tile. Therefore, if ^ is located anywhere but the top right corner, it will provide 
a net cost of at least +1, but if it is in the top right corner, it provides a net benefit: —1 to the total cost. 
Similarly, for Q, Q. and g], we arrange rules that cause them to be preferentially in the top left, bottom 
left, and bottom right corners, each providing —1 to the total cost when in the correct positions, and at least 
+1 elsewhere. The pairing of (Q, in either the vertical or horizontal directions has a weight 0. In the 
main layer of tiles, we use only the weights and +2, with +2 indicating a pairing that would be forbidden 
in the original TILING problem and indicating a pairing that would be allowed. Finally, we insist that 9, 
B' B' '^^ 9 special layer must correspond to the usual corner tile ti in the main layer. Then in 

order for there to be a solution with total cost at most —4, the new layer must in fact consist of Q, 
and ^ at the corners, and Q elsewhere in the new layer (as in figure |7]l, and the main layer is constrained 
in exactly the way it would be in the standard TILING problem. 

WEIGHTED TILING with periodic boundary conditions: We can of course apply Theorem 15. 2[ but 
when we allow weights, there is a simpler solution which avoids the caveats about the reduction. 

We now consider only odd N, not necessarily prime N. We will use the weights 0, +1, and +3, and set 
c = +2. Thus, to have a good enough tiling, we can have two pairings with weight +1, and all the others 
must have weight 0. 

As above, we add a new layer, again with five types of tile, B, and We assign a 

weight to any pairing with (Q, B) adjacent, either vertically or horizontally and with either orientation. 
However, we assign a weight +3 to have (Q, or d, B) adjacent vertically or horizontally. Thus, the 
preferred arrangement is for Q and ^ to tile the plane with a checkerboard pattern. Of course, they cannot 
do that completely when is odd. 

We allow Q to be above or below Q or ^ (weight for such pairings), and we allow to be to the 
left or right of Q or ^ (again, assign weight to any such pairing). To the right or left of we have 
weight to have another but above or below we have a cost of +3 to have or Similarly, 

above or below we have weight to have but we have weight +3 to have or ^ to the right 

or left of Thus, Q will preferentially form horizontal lines, and will preferentially form vertical 
lines. 

Finally, we have weight to have ^ above or below and weight +1 to have ^ to the right or left 
of 0. We have a weight +3 to have ^ to be adjacent to Q or ^ in any direction. Therefore, the only 
place ^ can be to provide total cost +2 is at the intersection of a horizontal line of Q's and a vertical line 
of ^'s, and there is no place it can go to have a cost of or +1. The cost of having a line of Q's and ITI's 
intersect without a ^ is even higher, as is the cost of not having a line of Q's or ^'s at all. The upshot 
is that any tiling of the new layer must be mostly a checkerboard of Q and B, with exactly one horizontal 
line of 0's and one vertical line of ^'s, which intersect in a Those horizontal and vertical lines will 
determine the boundary of an {N — 1) x [N — 1) grid, and we can set the boundary conditions at the corners 
via adjacency to the Q and ^ tiles. The arrangement is much Uke that of figure [6^, but with Q replaced 
by B and Q replaced by ■ 
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5.3 One-Dimensional Tiling 

In one dimension, the tiling problem becomes the following: 
Definition 5.5 1-DIM TILING 

Problem Parameters: A set of tiles T = {ti, . . . ,tm}- A set of constraints H C T x T such that if ti is 
placed to the left oftj, then it must be the case that {ti,tj) G H. A designated tile ti that must be placed at 
the ends of the line. 

Problem Input: Integer N, specified in binary. 

Output: Determine whether there is a valid tiling of a line of length N. 

We can also define a WEIGHTED 1-DIM TILING problem analogously to the WEIGHTED TILING 
problem in section [5^ 

Theorem 5.6 l-DIM TILING and WEIGHTED 1-DIM TILING are in P. 

Proof: Unweighted case: Let us create a directed graph with m nodes. The ith node corresponds to the ith 
tile type ti, and there is an edge between i and j iff {ti, tj) G H. We wish to know whether there is a path 
of length exactly N starting and ending at ti. Certainly this can be done in time poly(A/^), but we actually 
wish to do it in time polylog(A^). We will also generalize slightly to allow the left end of the line to have 
tile to and the right end to have tile ti. 

First, let us make a table of all paths which start and ending at the same node and contain no additional 
cycles. Such a path is called a simple cycle. A path is a simple path if it is a simple cycle or it contains no 
cycles. Let M be the set of all simple cycles, and let l{p) be the length of the path p. Clearly when p is a 
simple path, l{p) < m + 1, which is constant, so constructing the table of simple cycles takes constant time. 

Given any path p, we can decompose it into a simple path and a multiset P of simple cycles by progres- 
sively removing simple cycles from p until we are left with a simple path, adding each removed cycle to 
P. 

Definition 5.7 A multiset P of simple cycles is allowed /or the simple path p if there exists a path q which 
can be decomposed as above into exactly the multiset P plus p. 

Claim 5.8 The multiset P is allowed iff the underlying set Qfor P is allowed. 

ProofProof of claim: If we have any path through the node ti, we can insert a simple cycle for ti to get 
another possible path. If g is a path which can be decomposed to the set Q, then every node on every simple 
cycle in P appears in q, so we can insert as many copies as we need of the simple cycles in P. 

Conversely, let g be a path for the allowed multiset P. We can express qasa. tree, with the root being the 
simple path p and each node being a simple cycle. Cycle c' is a child of c if they share a node and c' is not 
the parent of c. The tree can be obtained by iteratively taking two consecutive occurrences of a node ti along 
a cycle such that there are no repetitions between the occurrences of t , and identifying the two occurrences 
of V. We wish to modify the tree so that the resulting path q' has the same multiset P as its decomposition, 
but with the property that there is a subtree rooted at p containing just one of every simple cycle in Q. Then 
we can take the path q corresponding to that subtree, and q has the set Q for its decomposition. 

That is, we need that within the tree, for every r G Q, there exists an occurrence c of r such that the path 
from p to c in the tree does not pass through two occurrences of any simple cycle. Since we can move any 
node corresponding to a simple cycle for ti to be a child of any simple cycle that contains ti, we can put the 
tree in the desired format by moving one instance of each r £ Q \xp the tree until it lies off of the first cycle 
which contains ti. ■ 

Claim 5.9 There is a tiling of a line of length N from to to ti iff 3 simple path p, f : M ^ U {0} such 
that: 
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1. p = {ti,^,ti-^,. . . ,ti^), with io = 0, = 1, l{p) = L + l, 
2- L+l + j:^^,,fiq)l{q) = N, 
3. f^^{Z^) is allowed for p. 

Given an allowed multiset P of simple cycles, the function f{q) represents the number of times the 
particular simple cycle q appears in the multiset, which means that is the underlying set for P. To 

see what Claim [S!9] is true, observe that the resulting path contains N nodes, as desired. Conversely, given p 
and /, we can create an allowed multiset P for p by including each simple cycle q a number of times equal 
to f{q). 

There are only a constant number of simple paths p and allowed sets P of simple cycles. By running 
over {p, P) and applying claim 15.91 the 1-DIM TILING problem reduces to determining if there exists /' : 
P ^ Z+ U {0} such that EgeP f'{q)l{q) = n', n' = n - {L + 1) - E.^p l{q). (We take f'{q) = f{q) - 1, 
with / given by the claim.) With fixed {p,P), we must therefore solve the following problem: Given a 
set of m' positive integers (the set of distinct l{q) for q G P), do there exist non-negative integers bk 
(the corresponding f'{q), summed over any cycles with equal lengths) such that "^j^akbk = N'l This is 
a special case of the Unbounded Knapsack Problem [|Lue75ll where the set of objects is fixed and only the 
total cost allowed varies. 

Note that, if we allowed 6^ to be negative, this would just be answered by determining if N' is a multiple 
of gcd(ai, . . . , flm')- Indeed, if A'^' is greater than m' g, where g is the least common multiple of the a^, this 
is all we need to determine: If gcd(ai, . . . , am')\'r, then we can write r = flfeCfc, with the property that 
afcCfc > —g for all k. Thus, if N' = m' g + r, then N' = 'Ylk'^k{ck + g/o-k), providing our solution. 
Conversely, if gcd(ai, . . . , a^') then it is not possible that ^^i. akbk = N'. For N' < m'g, matters are 
more complicated, but since g < m! is a constant, we can answer those cases via a look-up table. 

Weighted case: The 1-DIM TILING problem remains easy even if we take a weighted variant. In this 
case, the directed graph we get is a complete graph, but the edges have a weight associated with them, and 
we want to determine the minimal cost path of length from to to ti. Now we need to count both the length 
and the cost of each simple cycle. That is, we need to learn if 3 simple path p, f : M ^ Z+ U {0} such that: 

1. p = {ti^^,ti^,. . ■,tij, with io = 0,iL = 1, l{p) = L + l, 

2. L+l + Y.,^,,f{q)l{q) = N, 

3. c{p) + Y.qeMf{qyiq) ^ c, 

4. is allowed for 

where c{q) is the cost of path q. 

Once we fix (p, P), we again get a set ak of positive integer lengths l{q) for q ^ P, but we also have 
Cfe, the costs c{q) of g G P. We can assume f'{q) = for all but one cycle g € P for each length — we 
should pick the one with the lowest cost c{q). Now we wish to find non-negative integers bk that minimize 
J2k ^kCk subject to J2k ^kbk = -^'> which is again a special case of the Unbounded Knapsack Problem. 
We now have simple cycles of every length up to m, so m' = m, and we again set g < ml to be the least 
common multiple of the (i.e., the 1cm of 1, ... , m). We create a look-up table to give an optimal path for 
N < mg. When is a multiple of g, note that we can find an optimal path by simply calculating which 
simple cycle has the lowest ratio Cfc/a^; then we just use that cycle g/a^^ times. 

For other values of N' = dg + r, with d > m, we just look at the optimal path for (m — l)g + r and add 
(d — m + l)g/ak copies of the optimal simple cycle. This will be an optimal path for N'. To show this, note 
that we can assume without loss of generality that an optimal solution will not have 5^ > g/aj^ for more 
than one value of k: Suppose there were two values ki, k2 for which bk > g/ak, with Cfe^/ofe^ < Cfcj/o/tj- 
Then we can shift bk^ 6^1 + 5/«fci and bk^ '-^ ftfca ~ 9/'^k2^ leaving the total length the same without 
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increasing tlie cost. Furthermore, using the same logic, we know that the one value ki for which bk > g/ak 
must have a minimal value of Ck/au- For N' > mg, there must actually exist such a value of ki, so the 
optimal path must have been formed in the way we have described. 



5.4 Additional Symmetry 

We consider two additional kinds of symmetry. If we have reflection symmetry, then if {ti,tj) € H, then 
{tj,ti) £ H as well, and if {ti,tj) G V, then {ti,tj) G V also. That is, the tiling constraints to the left and 
right are the same, as are the constraints above and below. However, if we only have reflection symmetry, 
there can still be a difference between the horizontal and vertical directions. If we have rotation symmetry, 
we have reflection symmetry and also (t j ,tj) G H iff (t j ,tj) G V. Now the direction does not matter either. 

Note that these types of reflection and rotation symmetries assume that we can reflect or rotate the tiling 
rules without simultaneously reflecting or rotating the tiles. For instance, if a tile ti has a pattern on it (such 
as 9) that looks different when it is turned upside down, then when we reflect vertically, we also could 
reflect the tile, producing a new tile R{ti) in this case). We could define a reflection symmetry for this 
type of tile too: {ti,tj) G H iff {R{tj), R{ti)) G H, etc., but we just get the same complexity classes as 
for the case with no additional reflection symmetry. This is because we can add an extra layer of tiles with 
arrows on them and put on a constraint that any adjacent arrow tiles must point the same direction. While 
either direction will work (or any of the four directions in the case of rotation symmetry), one direction ends 
up preferred in any given potential tiling, so by looking at the arrow tile in a given spot, we can effectively 
reproduce rules that have no reflection symmetry. Thus, we address the case where the rules have reflection 
or rotation symmetry without simultaneously reflecting or rotating the tiles. 

Theorem 5.10 When the constraints for TILING have reflection symmetry, there exist N^., No G U {oo} 
such that for even N > N^ or odd N > No, a valid tiling exists, while for even N < Nf. and odd N < No, 
there is no tiling (except for N = 1, when there is a trivial tiling). 

• When we have open boundary conditions, either N^ = No = oo or N^ = No = 1. 

• When we have the four corners boundary condition, either No = 3 or N^ = No = oo. 

• When we have periodic boundary conditions, either N^ = 2 or Ne = No = oo. 

• When we have rotation symmetry, N^ and No are computable. 

When we have reflection symmetry but not rotation symmetry, we have been unable to determine so far 
whether N^ and No are computable for the four comers and periodic boundary conditions, respectively. 
Note that we do not consider just one corner fixed for the boundary condition, since that would break the 
reflection symmetry. 

Proof: To prove the theorem, simply note that given a valid tiling of an x grid, with iV > 4, we can 
extend it to a valid tiling of an [N + 2) x [N + 2) grid. The main observation is that we can repeat existing 
patterns when we have reflection symmetry, because if AB is a legal configuration, so is ABAB. 

In order to extend a tiling, we can do the following: strip off the leftmost column and bottommost row, 
and replace them with duplicates of the next two rows and columns. We can fill in the corner by duplicating 
the bottom left 2x2 square once we have stripped off the rows. Then the original leftmost column and bot- 
tommost row became the new leftmost and bottommost column and row, with some duplication to lengthen 
them to the right size. (See figured) 

This strategy handles even more general boundary conditions than the three main cases we consider. In 
particular, it also works if the tiling rules on the sides of the grid are completely different from the tiling rules 
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Figure 8: Extending a tiling of a 6 x 6 grid to 8 x 8 when the tiling conditions have a reflection symmetry. 
The spaces on the right mark the repeated rows and columns. 




Figure 9: Extending the tiling of a corner to a tiling of a 3 x 3 grid with fixed corner boundary conditions. 

in the interior of the grid. When the tiling rules are the same on the sides as in the center, except perhaps for 
the comers, we can copy the outermost two rows and columns, so this strategy works for N > 2. 

Four corners boundary conditions: When N is odd and we have the four-comers boundary condition, 
we can look to see if there is a tiling of a 2 x 2 grid with the corner tile in the upper right corner. If so, 
we can duplicate the right column to add a column on the left, and then duplicate the top row to add a row 
on the bottom. (See figure |9]) This gives us a tiling of the 3x3 grid with all four comers correctly tiled. 
Conversely, if there is no 2 x 2 tiling containing one of the corner tiles, then there cannot be a tiling of any 
N X N grid with > 1, which means that No = Ng = oo. 

Open boundary conditions: In the case of completely open boundary conditions, we need only check 
a 2 X 2 grid to see if there is a valid tiling. If not, there cannot be a tiling of any size N > 2 either. If there 
is, we can extend it as in figure |9] to get a tiling of the 3x3 grid as well, and extend as in figure [8] to get a 
tiling of any size grid. 

Periodic boundary conditions: When we have periodic boundary conditions, a valid tiling of a 2 x 2 
grid with open boundary conditions gives us a valid periodic tihng of a 2 x 2 grid as well, so we can extend 
it to all even N. However, if we try to apply the strategy of figure |9] to extend it to odd N, we potentially 
ruin the periodicity, so we do not know if No is computable in this case. If there is no tiling of a 2 x 2 grid 
with open boundary conditions, then there can not be a tiling for any A'^ with periodic boundary conditions 
and No = Ne = oo. 

Rotation symmetry: When we have rotation symmetry (with any boundary conditions), it will suffice 
to compute N(, and No as the minimum even and odd lengths that allow us to tile a single side of the square. 
If a tiling of one side exists, we can use this same tiling on all four sides and then fill in the center using 
diagonal stripes of identical tiles, as in figure [TO] In the case where there are special rules for the boundary, 
we will need to tile a side plus the adjacent row/column as preparation, but this presents little additional 
difficulty. 

To tile a single side, we can use an approach similar to the previous 1-dimensional case. However, now 
matters are much simpler, since the graph is now undirected. There are thus always many size 2 cycles, so 
we need only find the minimal even- and odd-length cycles for ti. That sets an upper bound on N^, and No- 
It might be one of these can be made smaller, but that is straightforward to check as well. ■ 

The weighted case is more difficult, and currently, when there is only reflection symmetry, we only have 
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Figure 10: A tiling of a 6 x 6 square with rotation symmetry. 

a result for the case of periodic boundary conditions. We will show that something similar to theorem [5. 101 
holds in this case as well. 

Theorem 5.11 For WEIGHTED TILING with additional symmetry, we have the following results: 

• With reflection symmetry and periodic boundary conditions: There exist N^, Nq G U {oo} such 
that for even N > Ne or odd N > No, a valid tiling exists, while for even N, 2c < N < Ne and odd 
N, 2c < N < No, there is no tiling. N^ is computable. 

• With rotation symmetry and open boundary conditions: Either there exists computable Nq € such 
that tiling is possible for N > Nq and impossible for N < Nq or there exists computable Nq € 
such that tiling is possible for N < Nq and impossible for N > Nq. (Nq depends on the weights and 
maximum allowed cost c.) 

• With rotation symmetry and periodic boundary conditions: One of the following three cases holds: 
There exists computable Nq £ Z+ such that tiling is possible for all N > Nq, there exists computable 
Nq € Z"*" such that tiling is impossible for all N > Nq, or tiling is possible for all even N and there 
exists computable Nq G Z+ U {oo} such that for odd N > Nq, tiling is possible. 

• With rotation symmetry and the four-comers boundary condition: There exist computable N(.,No € 

U {cxd} such that either there exists a tiling for any even N > N^ or for no even N > N^, and 
either there exists a tiling for any odd N > No or for no odd N > Nq. 

Currently, we do not know if the No, Ng for the periodic boundary conditions with reflection symmetry are 
computable. 

Proof: Reflection symmetry and periodic boundary conditions: Let us suppose we have a tiling of 
the N X N grid with periodic boundary conditions with total cost c' < c. Let the rows of this tiUng be 
. . . , Rn, and because of the periodic boundary conditions, let Rn+i = Ri- Let w{Ri) = J2 w{ta, th), 
where the sum is taken over adjacent tiles ta and ti, in row i?„ and let w{Ri, Ri+i) = ^{ta, tf,), where 
the sum is now taken over adjacent tiles ta, h, where ta is in row Ri and th is in row Ri+i. Then c' = 
'}^i{w{Ri) + w{Ri, Ri-^-l)). Notice that if we duplicate two adjacent rows, say Ri and Ri+i, then the total 
cost becomes c' + w{Ri) + + 2w{Ri, Ri^i) (because of the reflection symmetry). Now, 

N N 

[w{Ri) + w{Ri+i) + 2w{R„ Ri+i)] = 2 ^ [w{Ri) + w{Ri,Ri+i)] = 2c'. (6) 

i=l i=l 

Thus, there must exist some i for which 

w{Ri) + w{R,+i) + 2w{Ri, R^+i) < 2c' /N. (7) 

When c' is a constant, for > 2c', that means there is some f for which ti;(i?j)+ti;(i?j+i)+2w(i?i, iij+i) < 
0, since the weights are integers. Then we can duplicate rows Ri and Ri+i without increasing the cost. 
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Similarly, there are two columns which we can duplicate without increasing the cost. Thus, we also have a 
tiling for the {N + 2) x (N + 2) square grid with cost c" < c' < c. The only difference from theorem [5 .101 
is that there might be some small exceptions with N < 2c which have valid tilings but cannot be extended. 
A'^e then is the smallest even value for N such that A'^e > 2c and there is a valid tiling of the A'^e x -^e grid of 
weight c' < c. Similarly, is the smallest value for N such that > 2c and there is a valid tiling of the 
No X No grid of weight c' < c. 

Note that this argument requires that the cost c be a constant. It fails if c = c{N) grows at least linearly 
with N. 

When N is even, we can easily compute the exact minimal cost achievable. Consider all possible 2x2 
squares of tiles, and compute the cost of each by adding the costs of the four adjacent pairs in the square. 
Given any possible tiling of the periodic N x N grid of total cost c', let Sa,b be the 2x2 square whose top 
left corner is in location (a, b), and let w{Safi) be the cost of Safi- Then ^ w{Sa,b) = 2c'. Therefore, the 
minimal possible cost achievable for the N x N grid is N'^w/2, where w is the minimal cost of any possible 
2x2 square. When N is even, this is actually achievable by repeating a minimal cost square A^/2 times in 
each direction. The squares with their top left corner in location (a, b), with a, b even are exactly the square 
that we chose to repeat, and the squares in other locations are reflections of that square, which have the same 
cost. 

Rotation symmetry and open boundary conditions: This case is essentially trivial. We find the 
minimal weight between any pair of tiles, and simply tile the square in a checkerboard pattern with those 
two tiles. That is certainly the minimum cost achievable. 

Rotation symmetry and periodic boundary conditions: Find the lowest weight of a cycle of length N 
for the one-dimensional WEIGHTED TILING problem. Call that w. The minimum achievable total cost for 
2-D WEIGHTED TILING with rotation symmetry and periodic boundary conditions is then 2Nw {N rows 
and N columns each of cost w). This can actually always be done using the same tiling as the unweighted 
rotation symmetry case (figure [TOl). Thus, the problem reduces to the one-dimensional case, which we have 
argued is in P. However, with the additional reflection symmetry, the argument can be simplified. If all 
the weights are positive, then w > N and tiling is only possible for small N. If there exists a negative 
pair of weights, then for sufficiently large N, the minimal cost path will use many of that negative pair, 
and tiling is always possible for sufficiently large N. If the smallest weight is 0, then any sufficiently long 
one-dimensional cycle with constant weight w must use a 0-cost pairing, and it is possible to duplicate the 
pair of tiles used in that 0-cost pairing. Thus, the minimum weight cycle of length + 2 is also at most w. 
Therefore, we need only determine the minimum odd A^ for which w is 0, which is the length of the smallest 
odd-length cycle involving only zero-cost pairings. (When A^ is even and there is a zero-cost pairing, a 
zero-cost cycle is always possible by alternating between the two tiles involved in the pairing.) 

Rotation symmetry and 4-corners boundary condition: This case is more difficult, but still compu- 
tationally easy (in P). For lai^ge A^, the idea is that we will again tile most of the square using a pair of tiles 
ti, tj with minimal cost to be adjacent. That determines the asymptotic scaling of the total cost for large A^. 
However, the details are more complicated because of the effect of the comers. 

First, note that if the minimum weight for any pair of tiles is positive, then it is certainly not possible 
to tile an A^ X A^ grid when 2N{N — 1) > c. Conversely, if there is a pair of tiles {ti,tj) such that the 
cost for having them adjacent is negative, then we can always tile an N x N grid for sufficiently large A^ 
by taking a checkerboard of ti and tj, with only the comers different (if the designated comer tile ti is not 
already ti or tj). Suppose that w{ti,ti),w{ti,tj) < w. Then "sufficiently large" A^ means A^ such that 
2N{N - 1) - 8 > 8u; - c. 

The only difficult case is thus when there are pairs of tiles {ti,tj) with w{ti, tj) = 0, but no pairs with 
a negative cost. Consider the graph whose nodes are tile types ti and edges for any pairing of tiles with 
cost. Let Za be a set of tiles comprising a connected component of this graph with more than one tile, and 
let Z = UaZa- We will show that any tiling of total cost at most c must consist mostly of tiles from one 
particular Za, with only a constant number of different tiles, and that indeed it is sufficient to place all the 
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Figure 11: The shaded tiles are an upper-left 4-square. The tiles in the darker shade are the upper-left 
4-border. 

non-Za tiles near the comers of the grid. Then by considering all possible low-cost combinations of tiles 
near the corners, we will be able to determine the minimum achievable cost for the grid. 

We will call the k-square for the upper left comer to be the k x k square of tiles in the upper left corner 
of the grid. The upper-left k-border consists of the rightmost column and bottom row of the /c-square for the 
upper left comer. We make similar definitions for the other corners: the upper-right /c -border is the leftmost 
column and bottom row of the /c-square for the upper right corner, the lower left /c -border is the rightmost 
column and top row of the /c-square for the lower left comer and the lower right /c -border is the leftmost 
column and top row of the /c-square for the lower right comer. An example is given in Figure [TT] We say 
that a square is valid if its border contains only zero-cost adjacent pairs and the corner of the square that is 
also the corner of the larger grid has the correct corner tile. Since a A;-border must be connected that means 
that the edges within the border are all contained within a particular connected component Z^. 

We will determine if there is an iV x tiling of the grid by considering all valid tilings for the squares 
in each comer of the grid where we allow the sizes of the squares to go up to 2c. For each such combination, 
we require that the total cost of the tiling be some c' < c and that the tiles in each border come from the same 
connected component Z^. We then must determine whether it is possible to fill the remainder of the grid 
with zero-cost adjacent pairs, possibly interrupted by a few additional adjacent pairs with total cost at most 
c — c' (although it will turn out that this is never necessary). Suppose we have such a tiling, and consider tile 
ti G Za on the border of one square, and tile tj G Z^ on the border of a different square. Then, as shown 
below, there must exist a path from ti to tj that only includes zero-cost adjacent pairs of tiles. Note that the 
parity of the length of the path does not depend on which path we choose. 

More generally, given two /c-squares Ii and I2 for different comers of the grid, we are interested in 
whether it is possible to create a path of zero-cost adjacent pairs to get from ti e Ii to tj G I2, without 
necessarily being concerned about whether we can tile the whole grid: 

Definition 5.12 Suppose we have two k-squares Ii and I2 that are located at different corners of the grid, 
with borders from Z^- We say they are compatible for and Z^ if given tile ti on the border of Ii and tile 
tj on the border of I 2, there is a path of zero-cost pairings from ti to tj of the correct length for Ii and I2 to 
be placed together on an N x N grid. We say a set of four squares Ii, I2, I3, and 1^ which fit into the four 
comers of the grid is compatible for N and Z^ if they are compatible in pairs. 

Note that it does not matter which tile we pick from the border to define compatibility. This is because if 
we use t[ on the border of Ii instead of ti, we can take a path from t[ to ti using only the zero-cost pairings 
in the border of Ii, and concatenate it with the path from ti to tj. Similarly, if Ii, I2, and /a are vaUd squares 
located at three different corners, then if Ii is compatible with I2 and I2 is compatible with Is, it follows 
that Ii is compatible with I^. Also, if Ji and I2 are compatible for A^, they are also compatible for + 2: 
Because of the reflection symmetry, we can repeat a pair of adjacent tiles in a path to lengthen it by 2. 

For large A'^, if there is a tiUng of the N x N grid with total cost at most c, then there must exist, for 
some a, valid squares in the four comers Ii, I2, 13, and I4 which are compatible for A^ and Za such that the 
total cost of Ii, I2, 13, and I4 is at most c and the sizes of the squares are between c and 2c. To see why this 
is true, note that there can be only c non-zero-cost adjacent pairs in the entire tiling. Therefore, as k ranges 
from c to 2c, it must be the case that for each corner at least one /c-border contains only zero-cost adjacent 
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Figure 12: Extending a square to make a corner triangle. The black tiles have nonzero cost to be adjacent to 
any other tile. All other pairings that occur in the figure are zero cost. 

pairs. Furthermore, for each such pair of borders, there are at least c + 1 disjoint paths from a tile in one 
border to a tile in the other border. This follows from the fact that the borders themselves have more than c 
tiles. At least one of these paths must contain only zero-weight pairings. Thus we know that for sufficiently 
large grids, if a there is a tiling of total weight at most c, then there must be four valid squares for the four 
corners which are compatible for N and Za- The converse is also true: 

Lemma 5.13 Let Ii, I2, I3, and I 4^ be four squares which are compatible for N and Z^, N sufficiently 
large. Then there is a tiling of the N x N grid using Ii, I2, I3, and I4 in the four comers, with only tiles 
from Za used in the rest of the grid. 

Proof of lemma: Assume without loss of generality that Ii is in the upper left corner, I2 is in the upper 
right corner, I3 is in the lower left corner, and I4 is in the lower right corner of the grid. 

We will extend each square to create an isosceles right triangle whose outer diagonal border is composed 
of just one kind of tile. This can be done by progressively duplicating the outer sides of the square, each time 
shifting by one space towards the edge of the grid, as in figure [12] Finally, we can add additional diagonal 
layers to make sure the outer diagonal border is any particular tile from Za, and to increase the size of the 
triangle as much as desired. 

We can define compatibility of the newly created comer triangles in the same way as for the squares. 
When squares / and J for two different comers are compatible for N and Za, then their corresponding 
triangles remain compatible when extended as described above. In particular, there is a path of zero-cost 
tile pairs that will let us tile the top row of the grid between Ii and I2. Then we can copy the tiling of the 
interval between triangles diagonally down and to the left, much as in figure [TO] describing the unweighted 
case. This fills in the upper left corner of the grid, up to the diagonal line defined by the leftmost tile of 
the triangle for I2 and the topmost tile of the triangle for 13. We can do the same thing for the lower right 
corner of the grid, tiling the right edge between I2 and I4, and copying this tiling diagonally down and to the 
left. This leaves only a diagonal strip between I2 and I^, and we can tile that using a checkerboard pattern 
consisting of the outermost tile ti for the triangles and any other tile from Za which has a zero-cost pairing 
with ti. See figure [13] for an example of a complete grid tiled in this way. 

■ 

Thus, we can use the following algorithm to determine what the lowest achievable cost is for large even 
and odd grid sizes, when the lowest-cost pairings are cost 0: List the low-cost valid squares for each corner 
up to size 2c and determine which sets of four are compatible for large even N and for large odd N. The 
lowest total cost for each case tells us the minimum achievable cost. 



6 Conclusion 

We have shown that a class of classical tiling problems and 1 -dimensional quantum Hamiltonian problems 
can be proven hard, even when the rules are translationally-invariant and the only input is the size of the 
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Figure 13: A16 x 16 grid tiled by extending triangles on the four comers. 





problem. While this result was motivated by the desire to see if it could be hard to find the ground state in 
some physically interesting system, it is true that the tihng problem and Hamiltonian problem for which we 
prove hardness are not themselves particularly natural. Still, given that very simple cellular automata can be 
universal, it seems quite possible that even some very simple tiling and Hamiltonian problems are complete 
for NEXP and QMAexp respectively. 

We have considered a number of variants of the classical 2-dimensional tihng problems, but there are 
a few cases in which we have not been able to determine the difficulty. In particular, when we have addi- 
tional reflection symmetry (but not rotation symmetry), weighted tihng remains open with open boundary 
conditions or boundary conditions at the corners. In a number of the other cases with reflection or rotation 
symmetry, we know the problem is easy, but do not know whether the algorithm we have suggested involves 
an uncomputable parameter. Another very interesting question is the difficulty of the quantum 1-DIM TIH 
problem with reflection symmetry, or for that matter, the quantum 2-DIM TIH problem with reflection or 
rotation symmetry. We have no results for those cases. 
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